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Summary. — Generalized parton distributions have been introduced in recent 
years as a suitable theoretical tool to study the structure of the nucleon. Unifying 
the concepts of parton distributions and hadronic form factors, they provide a com- 
prehensive framework for describing the quark and gluon structure of the nucleon. 
In this review their formal properties and modeling are discussed, summarizing the 
most recent developments in the phenomenological description of these functions. 
The status of available data is also presented. 

PACS 12.38.-t - Quantum chromodynamics. 
PACS 12.39.-x - Phenomenological quark models. 

PACS 13.60.-r - Photon and charged-lepton interactions with hadrons. 



1. Introduction 

The composite nature of the proton was made manifest by the discovery of its anoma- 
lous magnetic moment [1] and confirmed by the observation of the electromagnetic form 
factors by Hofstadter and coworkers [2] . Nowadays it is firmly believed that the internal 
dynamics of hadrons such as the proton and neutron (collectively indicated as nucleons), 
is determined by the strong interactions between quarks exchanging gluons, as governed 
by quantum chromodynamics (QCD). However, a detailed description of the nucleon 
structure is still missing because QCD can only be solved in the perturbative regime of 
short distance phenomena probed in hard collisions, whereas the soft part of the inter- 
action corresponding to the long-distance behaviour requires a nonperturbative and/or 
numerical treatment like, e.g., in lattice simulations. 

In order to probe the internal structure of the nucleon large energy (j/) and momen- 
tum (q) transfers are necessary as in the so-called deep inelastic scattering (DIS) and 



ultimately in the Bjorken regime, i.e. when both i/ and Q 



q — V become very 



large with fixed xb — /2p-q (or xb — /2Mmv for a nucleon with mass at rest). 
Under these conditions scaling occurs, i.e. structure functions parametrizing the inclusive 
DIS cross section become independent of at fixed xb- In the parton model [3] this 
phenomenon is interpreted as the incoherent elastic scattering off the partons with the 
Bjorken xb being just the fractional (light-cone, longitudinal) momentum of the struck 
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parton. At the parton level one may distinguish three kinds of quark distributions, i.e. 
the quark density q{xB), the helicity distribution Aq{xB), and the transversity AxqixB)- 
The first two arc well- known quantities: q{xB) is the probability of finding a quark with 
a fraction of the longitudinal momentum of the parent (fast-moving) nucleon, regard- 
less of its spin orientation; Aq{xB) gives the net helicity of a quark in a longitudinally 
polarized nucleon, i.e. it is the number density of quarks with positive helicity minus the 
number density of quarks with negative helicity, assuming the parent nucleon to have 
positive helicity; in a transversely polarized nucleon, the transversity At^Cs^b) is the 
number density of quarks with polarization parallel to that of the nucleon minus the 
number density of quarks with antiparallel polarization. Information on the last distri- 
bution is missing on the experimental side because A.Tq{xB) decouples from inclusive 
DIS and therefore cannot be measured in such a traditional source of information. Other 
processes, such as polarized Drell-Yan dilepton production, are better suited for accessing 
transversity. 

Experiments show scaling violation with a dependence due to the contribution of 
gluons taking part to the scattering process as active particles or being radiated by the 
initial and scattered quarks as described by QCD with the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations [4]. Thus, QCD fits to DIS data also 
determine the gluon distribution. A large amount of data over the years have provided 
us with a fair description of quark and gluon distributions (see, e.g., [5]). 

More generally, in the Bjorken regime one probes space-time correlations along the 
light-cone. Whereas inclusive DIS involves diagonal matrix elements of certain operators, 
thus allowing a probability interpretation in terms of distributions, a full knowledge of 
the correlations can only be achieved by considering also the nondiagonal matrix elements 
of the same operators. This is possible in exclusive processes under suitable conditions 
where one can factorize short- and long-distance contributions to the reaction mechanism. 
These nondiagonal matrix elements can be parametrized in terms of generalized parton 
distributions (GPDs). GPDs have been introduced in the past in different contexts (see, 
e.g., [6, 7]), but have raised a large interest in the hadron commimity only when their 
importance was stressed in studies of deeply virtual Compton scattering (DVCS) [8, 9, 10] 
and hard meson production [11] in connection with the possibility of factorizing their 
contribution [12] and gaining information on the spin structure of the nucleon [8]. Being 
related to nondiagonal matrix elements, GPDs do not represent any longer a probability, 
but rather the interference between amplitudes describing different parton configurations 
of the nucleon so that they give access to momentum correlations of partons in the 
nucleon. 

The finite momentum transfer to the proton makes a second space-time structure of 

the process possible. Whereas in inclusive DIS the partonic subprocess is the scattering 
of a photon on a quark or antiquark, in the case of GPDs the virtual photon can also 
annihilate on a quark-antiquark pair with transverse separation of order 1/Q m the 
proton target. 

The momentum transfer can also have a transverse component. This provides a 
powerful tool to study hadron structure in three dimensions, because in addition to the 
information on the (longitudinal) behaviour in momentum space along the direction in 
which the nucleon is moving (as in the case of ordinary parton distributions), they also 
give insights on how partons are spatially distributed in the transverse plane [13] (as in 
the case of elastic form factors). 

The GPDs can also be viewed as the generating functions for the form factors of 
the twist-two operators governing the interaction mechanisms of hard processes in the 
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deep inelastic regime. These generalized form factors do not couple directly to any known 
fundamental interactions, but can be studied indirectly looking at moments of the GPDs. 

The most peculiar example arc the form factors of the energy momentum tensor, which 
allow to access the total and orbital angular momentum of the nucleon carried by quarks 
and gluons. 

In the last ten years GPDs have been investigated in great detail from the theoretical 
point of view, and many review papers already have summarized the progressing status 
of their understanding and modeling [14, 15, 16, 17, 18, 19]. 

Experimentally, measuring GPDs in exclusive processes is a big challenge requiring 
high luminosity and resolution. Therefore only very recently first dedicated experiments 
have been planned and performed. 

This review is an attempt to provide a general overview on the most recent develop- 
ments in the phenomenological description of the nucleon GPDs. After a summary of 
properties of the GPDs in sect. 2 and their physical content in sect. 3, the main lines 
of research in modeling GPDs for the nucleon are presented in sect. 4. Results obtained 
for quantities describing the nucleon structure are discussed in sect. 5, and the status 
of experimental investigation is presented in sect. 6. Conclusions are drawn in the final 
section. 

2. — Definition and properties of generalized parton distributions 

2T. Generalized form factors. - In quantum field theory one can construct currents, 
i.e. Dirac tensors as bilinear combinations of F matrices (F — 1, 75, 7'^, 7''75, cr'''^) and 
Dirac fields, ■tpTtp. Depending on the selected F matrix and its Lorentz properties one 
can define scalar, pseudoscalar, vector, axial and tensor currents. 

Matrix elements of the above currents on nucleon states with initial (final) momentum 
(p'^) and covariantly normalized as {p'\p) = 2p^{2tt)^6{p' — p), are usually expressed 
in terms of form factors. For example, with ip tpqi^) foi' quarks of flavour q, the matrix 
elements of quark (electroweak) vector and axial currents are decomposed as 



(2.1) {p'\Mow,mp) = u{p') 

(2.2) (p'|Vi,(0)7^75V'9(0)b) = U{P') 



u{p), 
u{p), 



where u{p) is the Dirac spinor normalized as u{p)u{p) = 2M]^, and = p'^ — and 
t = A"^. For each separate flavour q the vector current involves the contributions 
and F2 respectively to the Dirac and Pauli electromagnetic form factors of proton and 
neutron, F['" and i^f'", while the weak axial current involves the contributions g\ and 
gp respectively to the axial and induced pseudoscalar form factors. 

Restricting oneself to up and down quarks only and referring to their contribution in 
the proton, the form factors ^^e related to the physical 2* as 

(2.3) f'"2 = 2Ff^2 + J^i"2, Fl, = Fl, + 2Fl2. 

For the axial vector form factor one uses the isospin decomposition 
(2-4) 9A = yA + h9A, 9A = -h9A + y% 
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where gA {g\) is the isovector (isoscalar) axial form factor with 5a(0) = 1.267 and, 
from quark models, g\{t) = ^gA{t). Relations similar to (2.4) hold for the induced 
pscudoscalar form factor containing an important pion pole contribution through the 
partial conservation of the axial current. 

The above currents can be considered as particular cases of more general operators in 
QCD. In the deep inelastic regime of hard processes, where an operator product expansion 
is performed in order to overcome the problem of light-cone singularities arising as a 
consequence of the explored short distances, the major contribution comes from the so- 
called twist- two tensor operators (sec, e.g., [20]). Formally, the twist t is defined as the 
dimension d in mass units minus the Lorentz spin s of the operator, t = d — s. A spin-s 
tensor transforms as an irreducible representation of the Lorentz group. The maximal 
spin for a given number of Lorentz indices is achieved when they are all symmetrized. 
The irreducibility implies that the reduction to lower-spin tensors is not possible: as a 
consequence, the contraction of any pair of indices with the metric tensor gives zero. 
Thus, the Lorentz structure has to be traceless. 

In QCD, there are six towers of twist-two operators forming totally symmetric repre- 
sentations of the Lorentz group [14, 18, 19]: 



(2.5a) 








(2.56) 


q 






(2.5c) 






■■iV''-^Hq , 


(2.5d) 




= F^'^°'i V '^^ ■ ■ ■ 




(2.5e) 








(2.5/) 




= F^^^i V ■ ■ ■ 





where F"!^ = ^e°''^'^^Fjs is the dual field strength tensor with eoi23 = 1. In Eqs. (2.5) all 
indices within (• • •) are symmetrized and traceless, and 

(2.6) = i (p^ - 5^) , Vy, = ~di,-igtaA''^{z), T>i, = *dy.+igtaAl{z) 

are covariant derivatives expressed in terms of the gluon vector potential A'^{z). The 
label a = 1, ... ,8 is the octet colour label, and the ta are (one half of) the Gell-Mann 
matrices for the triplet representation of SU(3), with 

(2.7) [ta,h]=ifabctc- 
The gluon field tensor is 

(2.8) F;, = d^Al - d.Al + g UcA^Al, 

where g is a constant representing the coupling strength between ■^g and A^. In the 
following when not necessary the colour index will be omitted. 

Like the form factors of the electromagnetic current, additional information about 
nucleon structure can be found in the (generalized) form factors of the twist-two operators 
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when the matrix elements are taken between states of unequal momenta. Using Lorentz 
symmetry and parity and time reversal invariance, one can write down all possible form 
factors of the spin-n operator c>wi- -^^n-i -^^ ^2.5a) [14, 18] 



(2.9) 



n-1 



u{p') 



■ even 

n-1 . 

uip), 



A^'^A'^i • • • A''"-^) C«(i) Mod(n + 1, 2) 



where P = ^{p + p') is the average nucleon momentum, and Mod(n + 1,2) is 1 for even 
n and for odd n. Thus C^{t) is present only when n is even, and in general there are 
n + 1 form factors [21, 22]. 

Similar decompositions in terms of generalized form factors are possible also for nu- 
cleon matrix elements of the other twist-two operators in Eqs. (2.5). The problem of 
counting the number of allowed form factors has been addressed in Rcf. [21] making use 
of a method based on partial wave formalism and crossing symmetry, i.e. by considering 
the number of independent amplitudes in the crossed channel (pplOjO) corresponding to 
proton- antiproton creation from the twist-two source. In the case of the vector operator 
one finds exactly n + 1 form factors, as quoted above. The situation is more complicated 
in the case of the axial vector and tensor operators [22, 23]. 

As has been observed in Ref. [17], there is no C^(t)-like generalized form factor 
present for the axial vector, and the parametrization reads [22] 



(2.10) 



(p'\ ^^(0)7^5 i-D"' ■ ■ ■ iV^-'^g{0) \p) 

n-1 

= u{p') {7^^75 • • • A^-P'^'+i • • • P'^"-i)i«_,(i) 



i=0 
even 



+ 75 77i^A(^A''i • • • A'^-P^'+i • • • p/'-Os' ,(i) j u{p), 



for a total of 2 [2^] -|- 2 independent form factors. 



There is a total of 2[^^] + n + 2 independent form factors for the parametrization 
of the tensor operator ^ Eq. (2.5c) [22]: 



(2.11) {p'\ iPg{0)ia^'^''W' ■ • •ir>^''-^Vg(0) \p) 



= A S u{p') 

[ixv] (vni...) 



J2 iCT'^^'A"! • • • A^^P^^+i • • • Pi^^-' AT„,i{t) 



i=0 
even 
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n-1 



A^^i . . . A'^'P'^'+i • • • P^'"-^ATn,i{t) 



even 



n-1 




2Mjv 



i=0 



even 




n-1 



A^i 



. . . A^-P^-+i . . . P^— iBT„,i(t) 



u{p), 



odd 



where one has first to symmetrize and then to antisymmetrize as indicated [24]. 
Analogous form factor decompositions are possible for the gluon operators [17]. 

2'2. Generalized parton distributions. - In hard scattering processes, where hadrons 
and partons move fast in the ^-direction, it is natural to work with hght-cone coordinates 
in terms of two hght-hke four-vectors n+ = (1, 0, 0, 1) / \/2 and n_ = (1, 0, 0,-1)/ \/2, i.e. 



where u+ = w • n_ = + v'^)/\/2, v~ = v ■ n+ = (w° — w^)/-\/2 and v±_ = (0, v^_^Q). The 
relevant momenta are then the light-cone plus-momenta p+ and the relevant F structures 
become, e.g., rj,- = ^ ■ ri- =7+. 

Light-cone bilocal operators arising in hard scattering processes can be expanded in 
terms of the above currents making use of the relation 



and its analogs for operators involving other F matrices than 7+. In general, with 
the covariant derivative T) instead of d on the right-hand side, a Wilson link operator 
5^;, 5^;] appears between the operators at positions —^z and ^z, where 



and V denotes path ordering between a and h. In the light-cone gauge, = 0, the Wil- 
son link reduces to unity, as it will be assumed in the following. However, the condition 
j4+ = does not remove all gauge freedom because z~ -independent gauge transforma- 
tions are still possible with consequences on the correlation functions in hard processes 
sensitive to transverse momenta of partons [25, 26, 27, 28, 29]. Technical details needed 
in the calculation of the gauge link corresponding to a given partonic subprocess and a 
calculational scheme are provided in Rcfs. [30, 31]. 

Parton distributions are just defined in terms of matrix elements of light-cone bilocal 
operators between proton states of equal momenta [32]. In general, with initial (final) 



(2.12) 




2=0 



(2.14) 
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M + lA 



Fig. 1. (a) Kinematic variables in the symmetric frame of reference; (b) parametrization of 
the GPD H''{x,^,t) in terms of momentum fractions. 



momentum p {p') and helicity Ajv (A^) one defines a set of generalized quark distributions 
for a hadron with spin i that have been classified in Refs. [33, 34]: 



(2.15) 



(2.16) 



(2.17) 



J 47r 



vP+z 



2P+ 



2M, 



JV 



+=0, zj_=0 

u{p, Xn), 



^ e-^^^" (p', a;,! M-k^) 7+75 V,(|^) b, A;v) 



u{p',X'n) 



dz 



2P+ 

AxP+z~ 1^1 \l 



2M, 



AT 



2 + =0, Zj_=0 

Aiv), 



(p',A';vlV'g(-5-2)«C^+V9(|-2) \P,Xn) 



!+=0,2 



2P+ 



m(p', A^) 



P+A' - A+P' 



7+ A' - A+Y 7+P' - P+y 



JV 



AT 



u{p, Xn)- 



Because of Lorcntz invarianco the eight GPDs m , E^i , m , E'f , H^, H'^, can 
only depend on three kinematical variables (Fig. 1), i.e. the (average) quark longitudinal 
momentum fraction x = k~^/P~^, the invariant momentum square t and the skewness 
parameter ^ given by 



(2.18) 



P+-P'^ 
p+ + p'+ 



A+ 
"2P+' 



In addition, as in the case of parton distributions, there is an implicit scale dependence 
in the definition of GPDs corresponding to the renormalization scale /x^, i.e. the scale at 
which the QCD operators in Eqs. (2.15)-(2.17) are understood to be renormalized. 

The eight GPDs are all required to be real valued as a consequence of time reversal 
invariance, with support in the interval x,^ £ [—1, 1]. 

Similarly for gluons one has [34] 



PT / ^e^^^"^>',A^lP+'(-i^)P>(i^) \p,Xn) 



z+=0,zj_=0 
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2P+ 



u{p', A^) 



H9 7+ + E<^ 



2M, 



u{p,Xn), 



(2.20) I ^e-^^^"(p',A^|F+X-^^)i^.+ (|^) \pAn) 



z+=0,zx=0 



2P+ 



Ha ^+^5 + £;s 



2M, 



AT 



u{p,Xn), 



(2.21) 



P+ / ItT 



gixP+z-y^ A^l SF+\-lz) F+^ih) \p, A 



1 P+A^-A+P^ , 

S , , „, w(p , A;v) 



2P+ 2MArP+ 



2M, 



AT 



N I 

z+=0, Zi=0 

P+A* - A+P* 



„7+A'-A+-' ^ -+P'-P+7* 



M. 



N 



u{p, Xn), 



where a summation over i = 1, 2 is implied and S denotes symmetrization in i and j and 
subtraction of the trace. 

Quark and gluon operators can be usefully rearranged in order to make explicit their 
action on the parton helicity, as indicated in Tab. I. The corresponding matrix elements, 



(2.22a) Al,. 



-I 



dz 



z+=0,zj.=0 



1 



(2.226) A\, , , = T^T / e- 



dz 



ixP'^z /„/ \l 



(P'>AW|0^',^(^) \pAn) 



Z + =0,Z_L=0 

are similar to helicity amplitudes [17]. In particular they share the same property from 



Table I. - Quark and gluon operators, C^/^ and O^/^ respectively, for leading-twist GPDs and 
the parton helicity transitions they describe. It is implied that the first field is taken at —^z 
carrying helicity fx and the second at carrying helicity /u. 







X < 




X 


<ICI 


x>e, 






O", 

MM 
















iV;,7+(l + 75)V,, 








+ 


+ 


+ 




i^,7+(l-75)v^. 


+ 


+ 


+ 










-iiV;,(a+i-ia+2)V., 


+ 




+ 


+ 




+ 




iiV',(a+i+ia+2)V, 




+ 






+ 
















+ 


+ 


+ 




^[F+'Fi+ +iF+'Fi+] 


+ 


+ 


+ 












+ 




+ 


+ 




+ 




- + + iF+^F''+] 




+ 






+ 
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parity invariance, 

(2.23) A_v -^',-A„-M = (-1)^'--'''-^-+'' 



A 



Using light-cone helicity spinors explicit calculation [17] gives the matrix elements 



(2.24a) 

(2.246) 

(2.24c) 
(2.24rf) 



-+,-+ 



H + H ^2 ^ + ^ 
2 1 -^2 2 

H-H ^2 E-E^ 



1-^2 2 



++ 



VUT^ E-jE 
2Mn 2 ' 
VUT^ E + ^E 



2M. 



N 



for both quarks and gluons, where (p is the azimuthal angle of the vector D = — 
- 15/(1 + 0, i-c. eJ'fi = (D^ + iD^)/\D\. Therefore, the four GPDs H, E, H, E are 
partori helicity conserving and chiral even. The distributions H and E are sometimes 
referred to as unpolarized and H and E as polarized because H and E correspond to the 
sum over parton helicities, and H and E to the difference. 
In the case of parton-helicity flip, for quarks one obtains 



(2.25a) ^++,+_ 
(2.256) 

(2.25c) A%__ 
(2.25d) 



2M, 



2 



2M 



JV 



'9 

T ) ' 



^0 - t frq 

AMI ^' 



where 

(2.26) -h = 

is the minimal value for —t at given ^ . Analogous expressions hold for gluons with 
an additional global factor e"^ i/l — ^2 — tj {2Mn) on the right-hand side. Other 
helicity combinations are given by parity invariance. Therefore, the four GPDs Ht, Ht, 
Ej" , Ej" are parton helicity flipping and chiral odd. 

By inverting the set of Eqs. (2.24), (2.25) the different GPDs can be separately ex- 
tracted from the amplitudes ^A'jy/ii'.Aw/i- 



10 



SIGFRIDO BOFFI and BARBARA PASQUINI 



2'3. Forward limit. - In the forward case, p = p' , both A and ^ are zero. As ^ ^ 0, 
also X ^ xb, where the Bjorken variable xb is the fraction of the longitudinal momentum 
carried by the active parton. In this case the functions H'^, H'^ and reduce to the 
usual DIS parton distribution functions, i.e. 



(2.27) 
(2.28) 
(2.29) 



iJ«(xB,0,0) = 
#«(a;s,0,0) = 
jy|(:cB,0,0) = 



q{xB), Xb > 0, 

-q{~XB), XB < 0, 

Ag(xB), XB > 0, 

Ag(-XB), XB < 0, 

ArqixB), XB > 0, 
-Arqi-XB), Xb < 0, 



where q{xB), {q{xB)), Aq{xB) {Aq{xB)) and ATq{xB) {ATq{xB)) are quark (antiquark) 
density, helicity and transversity distributions, respectively: 



(2.30) q{xB) = / ^ e»-f^^" (^1^5(0) 7+ '4>{z)\p) 



2+=Z I =0 



(2.31) Aq{xB)= y'^e--f+^"(p5|| 1^(0) 7+75 ^-(-2)1^^11) 



(2.32) ATq{xB) = /^e"-^''^ (p5^|V'(0) 7V75 V'(^)|p5±> 



z+=«j.=0 



with 5*11 {Si_) being the longitudinal (transverse) nuclcon-spin projection. 

No corresponding relations exist for the functions i?"^, E'^ , and because in the 
forward limit they decouple in their defining equations. However, they do not vanish. 
In particular, i?'(.T,0,0) carries important information about the quark orbital angular 
momentum (see. Eq. (2.67)). In contrast, E^ vanishes identically being an odd function 
of ^ by time reversal symmetry [34] . 

For gluons one has 



(2.33) HS{xB,0,0)=XBg{xB), H\xb,0,Q) = xb Ag{xB), Xb > Q. 



All the gluon helicity-flip matrix elements go to zero in the forward case corresponding 
to the collinear limit {t = to in Eqs. (2.25)) and therefore for a spin-i target decouple 
from any observable for collinear scattering. 

2'4. Polynomiality . - Introducing a light-like vector n'^ , which is conjugate to in 
the sense that P ■ n = 1, and contracting both sides of Eq. (2.9) with n^n^^ ■ ■ ■n^^_^, 
one obtains 



(2.34) 



.n^.---n^,.-.{p'\0'i"-'\p) 

= u{p') i u{p) Hl{i, t) + n(p') " u(p) El{i, t), 
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where 

n-l 

(2.35) H^^iU) = E <»W (26' +Mod(n + 1,2) C«(t)(20", 

»=0 

even 

n-l 

(2.36) EU^, t)=J2 KM (20^ - Mod(n + 1, 2) C?,{tm)\ 



i = 

even 



Recalling the relation (2.13) the functions H^{^,t) and E^{^,t) are easily recognized to 
be the n-th Mellin moment of the GPDs and £'^(x,^, t), respectively, i.e. 

(2.37) j ^dxx^-^H\x,^,t) = Hl{^,t), 

(2.38) j'^dxx^-'E'^{x,^,t) = El{^,t). 

Eqs. (2.35) and (2.36) show that the n-th MeUin moments of the GPDs Hi{x,^,t) and 
E''{x,£,,t) are polynomial functions of ^ with highest power n. This is called the poly- 
nomiality condition. In QCD the polynomiality property of GPDs follows from their 
hermiticity and Lorentz, parity and time-reversal invariance [14]. 

An alternative notation is often used instead of (2.35) and (2.36) [16], i.e. 

(2.39) /' dx X--' H\x, 1 1) = 4"^ (i) + 4"^ {t)e + ...+ [ ^EI ^" . " 

^ ' J-i ^ ' ° ^' ^ ^' \/ii'^i(f)e-^fornodd, 

(2.40) dx X"- E^ix, ^, t) = ei-\t) + et\t) + . . . + ( ^^^f " , f " 7; 
7-1 [e;^i(i)^" ^ for n odd, 

where flavour indices are suppressed for brevity. For a spin-^ particle the coefficients 
in front of the highest power in ^ for even n are related to each other and arise from 
the so-called D-term D'i{z,t) with z = x/^ [35, 36], which has finite support only for 
\x\ < 1^1 , according to 

(2.41) /i^("\t) = -e^("^(t) = /' dz z^-^D\z,t). 



Similar polynomiality conditions can be derived for the other twist-two operators. For 
example, in Eq. (2.10) time reversal invariance constrains m(p')75"(p) come with odd 
powers of A". One then has the following sum rules [17] 

/I n-l 
dxx^-^H'i{x,^,t)^ E(2?)M^^,(t), 
'1 i=0 

even 

/I n-l 
dxx-'E'i{x,i,t)=Y,m'BlS)^ 
-1 ITo 
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where the highest power in ^ is n — 1 instead of n. 
For chiral-odd GPDs one has 



(2.44) 
(2.45) 

(2.46) 
(2.47) 



J dxx''-^Hl{x,^,t) 

j dxx''-^E^{x,i,t) 

j dxx''-^H^x,^,t) 

j dxx"-^E^{x,^,t) 



n-1 



i=0 
even 

»=0 

even 

n-1 

=0 
en 

n-1 

5](20^S^„,,(t). 



»=0 

even 



i=0 

odd 



The polynomiaUty property has been shown to lead to integral relations of the form 



(2.48) 



1 



dxF{x,^,t) 

-- I dxF{x,-,t) 
J-i <^ 



1 



— a- 



U!^ — X L0^ + X 
1 



1 



- I{oj, t) 



for (T = ±1 and any a; > 1, where F is one of the GPDs [37]. The only cases with nonzero 
T{lo, t) occur for unpolarized distributions (F = H'', E'^) and a = +1, where I(u;, t) takes 
contribution from the D-term, i.e. 



(2.49) 



± X{u>, t) = 2 



L 



dx 



D'i{x,t) 

U} — X 



with the sign + (— ) in the l.h.s. to be taken for H'^ {E'^). Thus polynomiality ensures the 
possibility of writing dispersion relations for the GPDs encoding causality and analyticity. 
In order to account for divergent contributions of both valence and sea quarks fixcd-i 
subtracted dispersion relations have been proposed [38] with the subtraction constant 
defined by the D-tevm (see also [39]). 

2 "5. Sum rules. - Relevant sum rules are obtained by forming the first Mellin moment 
of quark GPDs. For n = 1 in Eq. (2.34), from (2.37) and (2.38) one identifies and 

El (or /Iq^^ and eQ^^ in the polynomial expansions (2.39) and (2.40)) with the Dirac and 
Pauli form factors appearing in Eq. (2.1), i.e. 



(2.50) j dxHi{x,^,t)=F^{t), j dxE'>{x,^,t)=F^{t). 

As a consequence of Lorentz invariance, correctly the integrals only depend on t because 
integrating over x removes all reference to the particular light-cone direction with respect 
to which ^ is defined, so that the result must be independent of ^. At t = the quark 
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form factors and are normalized so as to yield the normalization 2 for the up 
quark distribution in the proton and 1 for the down quark distribution, i.e. -^"(0) = 2, 
Ff{0) = 1. The quark form factors and F2 arc normalized through the proton and 
neutron anomalous magnetic moments = 1.793 and k" = —1.913, i.e. F2{0) = /t" = 
2kP + k" = 1.673, Fi{0) = k'^ = kP + 2k" = -2.033. 
Similarly, one has 



(2.51) 



dxH''{x,^,t)=g\it) 



dxE\x,^,t)=gUt) 



where g\ and gp are the axial and induced pseudoscalar form factors of Eq. (2.2). The 
normalization of g\ and gj, a,tt = can be derived using the isospin decomposition (2.4) 

and the corresponding one for gp. 

There are also tensor form factors associated with the first moments of Ht, Ht and 
Et, whereas 



(2.52) 

In particular 



dxE^{x,^,t) = 0. 



(2.53) 



dx 



2H^{x,0,0) + E^{x,0,0) 



= K. 



'T' 



where nip describes how far and in which direction the average position of quarks with 
spin in the x-direction is shifted in the y-direction in an unpolarized nucleon [40] . Thus 

Kp governs the spin-flavour dipole moment in an unpolarized nucleon and plays a role 
similar to the anomalous magnetic moment k"^ for unpolarized quarks in a transversely 
polarized nucleon. 

The second Mellin moment of GPDs is relevant for the spin structure of the nucleon. 
The angular momentum operator in QCD contains quark and gluon spin and orbital 
angular momentum [41] with an important role of the orbital angular momentum in the 
formation of the total spin of the nucleon [42, 43]. It can be written in gauge-invariant 
form as J = ^ Jq + Jg, with 



(2.54) 
where 



(2.55) 



9' 



^9 



p 



are the gauge-invariant quark and gluon parts of the (Belinfante improved) symmetric 

energy-momentum (stress) tensor. They are twist-two operators of the type and 
O^^ discussed in section 2.1, and their matrix elements can be expanded in terms of 
generalized form factors: 



(2.56) W^imv)=u{jJ) 



_|_ pfi 



i(PMa-P + P-a^P)A, 
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u{p), 



where the form factor c{t) enters the quark and gluon parts with opposite signs in order to 
account for conservation of the total energy- momentum tensor. Alternatively, by means 
of the Gordon identity 

(2.57) 2MNu{p'h"u{p) = u{p'){ia'"'A^ + 2P«)u(p), 
Eq. (2.56) can be rewritten as [16, 44] 

(2.58) {p'\o^:',mp) = u{p' 



N 



2M 



N 



with 

(2.59a) 
(2.596) 
(2.59c) 



Aim = Mr(t), 



u{p) , 



From Eq. (2.56) one immediately has the following sum rules for the GPDs: 

(2.60) j dxx [H''{x,tt)+E''{x,tt)]=Al,{t) + Bl,{t), 

(2.61) J dx [Ha{x,^,t) + ES{x,^,t)]=Alo{t)+Blo{t), 

where the ^ dependence drops out eliminating the contribution of Cq^t), and the dif- 
ferent power of X in the integrals reflects the different forward limit of quark and gluon 
distributions. Extrapolating the sum rules to t = 0, the angular momentum carried by 
each parton species is found to be [8, 10] 

(2.62) (J3) = i [A%{Q) + Bl,{Q)\ , (J^) = \ [Al^Q) + i?,%(0)] . 

Since A^'q (0) give the momentum fractions of the nucleon carried by quarks and gluons, 



(2.63) 



/ 



dx 



X ^il«(a;,0,0)-hfl'»(a;,0,0) 



5^Al,o(0) + Af,o(0) = l, 



the total B{Q) = i?| q(0) + i?f q(0) vanishes for any composite system [45, 46]. Then, 
for the total angular momentum of a proton moving in the f-direction and polarized in 
the helicity eigenstate Ajv = +1 one recovers the value (J^) = ^q{Jq) + (Jg) = \- In 
addition, according to an extension of the equivalence principle of general relativity to 
describe the interaction of the nucleon with the external gravitational field one arrives 
to the interpretation of -B(O) as an anomalous gravitomagnetic moment being the analog 
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of the anomalous magnetic moment [47] . There is also evidence supporting the conjec- 
ture that the equivalence principle is valid separately for quarks and gluons resulting in 
exact cquipartition of momenta and angular momenta in the nucleon. The most precise 
numerical support comes from lattice calculations [48]. 

The sum rules (2.62), known as Ji's sum rules, are renormalization scale dependent. 
By solving the combined evolution equations for quark and gluon helicity and orbital 
angular momentum one finds in the asymptotic limit [49] 

where n/ is the number of flavours. Thus the partition of the nucleon spin between 
quarks and gluons follows the well-known partition of the nucleon momentum [50] when 
it is probed at infinitely small distance scale. 

It would be desirable to further split the Ji's sum rules (2.62) into spin and orbital 
angular momentum parts as [41] 

(2.65) i = 5^[iS'' + iJ+A5 + i„ 
where 

(2.66) ^^"/ j d^^ii^) 

are the quark spin and orbital angular momentum expressed in terms of the helicity distri- 
bution Aq{x) = H'^{x, 0, 0) and the orbital angular momentum density Lq{x). Analogous 
definitions hold for the gluon spin and orbital angular momentum, Ajf and Lg, respec- 
tively. However, the matrix elements of the gluon spin are not invariant under general 
gauge transformations [51]. One can only achieve a gauge-invariant separation of the 
quark orbital angular momentum [52]; 

(2.67) Lq{x) = ^x [H^ix, 0, 0) + E^ix, 0, 0)] - |fl"«(a;, 0, 0). 

When dealing with transverse total angular momentum, since there is no gluon spin 
contribution one can establish the transverse spin sum rule as follows [53]: 

(2.68) \ = \T.fdx l^Tq{x) + ^ (L^), 

where /S.Tq{x) = H!^{x,0,0) is the transvcrsity distribution and L±_ is the component 
of the angular momentum along the proton spin, polarized perpendicular to the proton 
direction of motion. 

The quark angular momentum can be further decomposed with respect to quarks of 
definite transvcrsity making use of the second Mellin moments o' ^T2 o ^^'^ ^T2 o 
the chiral-odd GPDs and E!^, respectively. If Jq^^^ is the angular momentum 

carried by quarks with transverse polarization in the -l-x-direction in a proton polarized 
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along X, the expectation value of the transverse asymmetry 6^Jg = Jg — Jg is given 
by [40] 



(2.69) 



{S'^J^) = i A|.2,o(0) + 2i^2,o(0) + S«2,o(0) 



2'6. Positivity. - Positivity bounds have extensively been studied in [54, 55, 56], 
where one can also find reference to previous work. They are all based on the positivity 
of the norm in the Hilbcrt space of states. There is a hierarchy of inequalities relating 
GPDs to the three quark distributions, i.e. the unpolarized, polarized and transversity 
distributions: i) strong inequalities where GPDs are bounded by combinations of all 
three distributions; ii) weaker inequalities without the transversity distribution; Hi) still 
weaker inequalities where GPDs are bounded only by the unpolarized distribution. 

In particular, for the region |^| < |x| one has [55] 



(2.70) 



1-e 



.E^x,U) 



< 



q{xi)q{x2) 
1-e ■ 



The function q{x) is the usual (forward) distribution function for unpolarized quarks of 
flavour q taken at values 



(2.71) 



X + ^ X - ^ 

= -r— -r, X2 = 



1 + C 1-^ 
The bound (2.70) is stronger than the following inequality derived in [57]: 



(2.72) 



e 



:Ei{x,tt) 



< 



iq{xi)q{x2) 



Prom (2.70) one also finds the bound for E'^ derived earlier in [57] up to a typo there: 

2Mn 



(2.73) 



\Ei{x,^,t)\ < 



^/q{x^)q{x2). 



Similarly, 



(2.74) 
and 



Ei{x,^,t) 



< y/q{xi)q{x2) 



(2.75) 



H''{x,^,t) 



< 



I -t q{xi)q{x2) 

to-t 1-e ' 
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3. Physical content of generalized parton distributions 

3'1. Parton interpretation. - The physical content of the bilocal operators entering 
the definition of quark GPDs becomes transparent in light-cone coordinates decomposing 
the Dirac field into the sum of the 'good' and 'bad' components corresponding to the 
independent and dependent degrees of freedom, i.e. tl'q,+ = -f+V'q a^nd tpq,- = P-i'q: 
with projection operators P± = ^7^7^- 

The 'good' fields have the following Fourier expansion in momentum space [58] 

(3.1) ^,,+ {z-,z^)= I ^^e(fc+) Y^{b,,ik+,k^)u+{k,i,)e-^'^^^-+^''--- 

tj. 

and similarly for V'g,+ . The spinors u+{k,fi) — P+u{k,ij) and v+{k,iJ,) = P+v{k,y) are 
the projections of the usual quark and antiquark spinors with helicity /z. The quark (an- 
tiquark) creation and annihilation operators {S) and h (d) obey the anticommutation 
relations 

(3.2) {b,^^'{k'+,k'^),bHk+,k^)} = {d^.^,{k'^.k'^).^l\^{k+,k^)} 

= 2k+{2wf 5(/e'+ - k+) S{k'^ - ki_) Sg,g 5^,^. 

Then, for example, the bilocal quark field operator Vig(— ii;)7+i/'g(|^) can be written as 

a density operator in terms of 'good' light-cone components, i.e. V^ipg _^_{—^z)iljq^^{^z), 
and the operator on the left-hand side of Eq. (2.15) becomes [57] 

(3.3) J ^e-^^^-^,(-i^)7+^(i^) 

X {s{'2xP+ - k'+ - k+)bl^,{k'+,k'^)bg^{k+,k^)ul{k',t,')u+{k,iJ,) 

+ 5{2xP+ + + k+) d,^> (fc'+, k'^)d\^{k+,k^) 4(fc', ^l')v+{k, 
+ 5{2xP+ + - k+) dg^^ (fc'+, k'^)bg^{k+,k^) vlik', fi')u+{k, m) 
+ 5{2xP+ - + k+)bl^,{k'+,k'^)dl^{k+,k^)ulik',t,')v+{k,tx)] . 

Which of the four terms in (3.3) contributes is determined by the positivity conditions 

fc+ > and k'~^ > for the parton momenta, together with momentum conservation, 
which imposes k~^ — k'^ = p+ — p'~^ = 2^P+. Assuming ^ > 0, the interval x G [—1, 1] 
falls into three regions according to whether |x| > ^ or jxj < ^: 

1. for X £ [(,, 1] both momentum fractions x + ^ and x — ^ arc positive. The corre- 
sponding GPD describes emission and reabsorption of a quark; 

2. for X G [— ^, £] one has x + ^ > and x — ^ < 0. The second momentum fraction can 
be interpreted as belonging to an antiquark with momentum fraction ^ — x emitted 
from the initial proton; 



18 



SIGFRIDO BOFFI and BARBARA PASQUINI 



(1 + 0^'^- 



it b 



d S 



{i+i)p+ 



{i + QP^ 



-{X + i)p+ 



Fig. 2. - Diagrams showing ^ symmetry of GPDs: a) a; > ^, b) — ^ < a; < ^, c) x < — ^. 



3. for X e [—1,—^] both x + ^ and a; — ^ are negative. The corresponding GPD 
describes emission and reabsorption of antiquarks with momentum fractions ^ — x 
and — ^ — X, respectively. 

The three situations are iUustrated by the diagrams on the left-hand side of Fig. 2. By 
assuming ^ < a similar analysis gives the diagrams on the right-hand side of the same 
figure [14]. 

Diagrams (a) and (c) in Fig. 2, arising from the b^b and d^d terms in Eq. (3.3), 
thus generalize the situation illustrated for the ordinary parton distributions. In these 
regions GPDs will evolve according to modified DGLAP equations [4]. Therefore these 
are called DGLAP regions. Diagram (b) in Fig. 2, corresponding to the middle region, 
— ^ < a; < ^, does not have a counterpart in parton distributions. This diagram arises 
from the db term in Eq. (3.3) and corresponds to the emission of a quark-antiquark pair 
with momentum — A. The middle region is thus similar to that in a meson amplitude and 
hence will evolve according to modified Efremov-Radyushkin-Brodsky-Lepage (ERBL) 
equations [59]. In this ERBL region, GPDs contain completely new information about 
the nucleon structure, because this region is completely absent in DIS, which corresponds 
to the limit ^ ^ 0. 

3'2. Impact parameter space. - GPDs are defined in terms of amplitudes in momentum 
representation involving proton states with a sharp plus-momentum. One can also use 
a mixed representation keeping momentum in the light-cone plus-direction and Fourier 
transforming from transverse momentum to transverse position. In such a representation 
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hadron states with definite plus-momentum can be localized in the transverse plane at 
a definite position to be interpreted as an impact parameter [60, 13]. This is possible 
according to the uncertainty principle and as a consequence of the fact that transverse 
boosts are purely kinematical, i.e. in the light-front framework they form a Galilei 
subgroup of the Poincare group. The position in the transverse plane coincides with the 
'center of momentum' R± = '^ikfb^ /X^j^i^ '^^ partons in the proton, given in 
terms of their plus- momenta kf and transverse positions (impact parameters) bj. 

In the limit ^ ^ 0, but t ^0, i.e. when the momentum transfer is purely transverse, 
a two-dimensional Fourier transform gives the distribution [61, 62] 

(3.4) W{x, b^) = [ ^ e-^"-^ H'^ix, 0, - A^), 



{2^r 

that depends on the impact parameter b only via its square thanks to rotation invariance. 
In addition, as a consequence of the Fourier transformation the matrix element is now 
diagonal in the plus-momentum and the impact parameter of the proton states that can 
be taken at R± = 0. Therefore, the impact-parameter dependent distribution (3.4) can 
be interpreted as the probability density to find a quark with momentum fraction x at 
transverse distance 6 from the center of momentum R± . 

Similar considerations can be done for the other GPDs. In particular, in order to 
develop a probabilistic interpretation for 

(3.5) S'^ix, b') ^ J f^e-'^-^E'^ix, 0, - A^), 

it is necessary to consider helicity flip amplitudes because otherwise E{x, 0, t) does not 
contribute. 

When A+ = for polarized nucleon states one has 

(3.6) J ^e-^^^" + A, t I Vi« (0) (^-) | P, t) = H'^{x, 0, - A^), 

(3.7) I ^e"^"^- {p + AA\i^AO)lUA^-)\pd)= -^^^^E^ixA-A'). 



Therefore, taking the nucleon polarized in the f-direction, one finds 

W) 



(3.8) p,(x,6)= / ^e-^-" 



ffn^,0,-A^) + ^i;^(x,0,-A^) 



Due to the contribution of S'^ in Eq. (3.8) that breaks the axial symmetry of 7i^, the 
quark distribution in a transversely polarized proton is distorted sideways in the trans- 
verse plane. This transverse distortion, combined with attractive final-state interactions, 
may give rise to the so-called Sivers effect [63] and relatively large transverse single-spin 
asymmetries [61, 64, 65, 66]. 

With probes exploring the nucleon on the scale of l/Q ^ 1 fm, the impact-parameter 
representation thus offers the possibility of performing a real femto-photography of the 
interior structure of the nucleon [67]. 
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This kind of analysis can be extended to the case of nonzero ^ [68]. In this case 
GPDs in impact parameter space probe partons at transverse position b with the initial 
and final state proton localized around R± = but shifted relative to each other by 
an amount of order ^b. In the DGLAP region the impact parameter gives the location 
where a quark or antiquark is pulled out of and put back into the proton. In the ERBL 
region the impact parameter describes the transverse location of a quark- antiquark pair 
in the initial proton. 

An extension of this framework to give a quark imaging in the proton has been pro- 
posed in Rcf. [69] developing the concept of the quantum phase-space Wigner distribution 
for quark and gluons in the rest frame of the proton and relating it to the transverse- 
momentum dependent parton distributions and GPDs. 

Otherwise, the Fourier transform of the DVCS amplitude with respect to the skcwncss 
parameter C = 2^/(1 -|-^) at fixed momentum transfer t provides an image of the nucleon 
in a boost-invariant variable a, the coordinate conjugate to light-front time [70]. The 
results obtained in a simple rclativistic model for spin-i systems are analogous to the 
diffractive scattering of a wave in optics where the distribution in a measures the physical 
size of the scattering center in a one-dimensional system. 

4. - Modeling GPDs 

It is quite difficult to calculate GPDs from first principles of QCD. Lattice simulations 
have received increasing attention in recent years and promising results have appeared. 
However, we are still far from explicitly calculating GPDs on the lattice. Deeply virtual 
Compton scattering and hard exclusive electroproduction of mesons give a theoretical 
possibility of experimentally constrain GPDs. However, this is a quite difficult task since 
observables involve convolution of the GPDs with hard scattering coefficients and not 
GPDs themselves. Therefore, there is continuing interest in modeling GPDs. 

There are basically two approaches. One is using ansatze to parametrize GPDs to 
be used in phenomcnological analyses. Here, the most popular choice is to parametrize 
the hadronic matrix elements which define GPDs in terms of double distributions [71, 
72], modeled by assuming a factorized t-dependence determined by some form factors. 
However, since this assumption is not strictly valid, alternatives have been suggested 
with the so-called dual representation [73], with constraints derived from data on the 
first Mellin moments, i.e. form factors [74, 75], or from simultaneous fits of data and 
lattice calculations of the higher Mellin moments [76]. 

Fitting procedures require suitable ansatze as well as theoretical developments. A 
promising approach along these lines has been discussed quite recently combining dis- 
persion relation and operator product expansion techniques and deriving the conformal 
partial wave decomposition of the DVCS amplitude in terms of complex conformal spin 
to twist-two accuracy [77] . The time-ordered product of the two electromagnetic currents 
sandwiched between the initial and final hadronic states in the Compton tensor is ex- 
panded in the basis of so-called conformal operators [78] . Based on the analyticity of the 
Compton tensor the conformal operator product expansion can be combined with the dis- 
persion relation technique to express the Taylor expansion of the Compton form factors 
with respect to u = 1/^ in terms of the Mellin-Barnes integral representation with coef- 
ficients that are independent of the renormalization/factorization scheme. With suitable 
ansatze for the conformal GPD moments in terms of hadronic partial wave amplitudes 
in the i-channel the Mellin-Barnes representation is adequate for building fiexible fitting 
procedures of DVCS and hard meson production [79] . 
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The other approach is a direct calculation using effective quark models. After the 
first calculation within the MIT bag model [80], GPDs were calculated in the chiral 
quark-soliton model [81, 82, 83, 84], the Nambu-Jona-Lasinio model [85], the light-front 
Hamiltonian approach [86], the nonrelativistic [87] and light-cone [88, 89, 90, 91] con- 
stituent quark model, the meson-cloud model [92] . 

After a brief account of the status of lattice simulations, in the following subsections 
the main lines of research within the two approaches are briefly reviewed. 

4'1. Lattice simulations. - Lattice QCD offers a unique opportunity to calculate Mellin 
moments of GDPs from first principles. The first investigations of GPDs including studies 
of the quark angular momentum contributions to the nucleon spin have been presented 
by the QCDSF collaboration in quenched QCD [93] and by LHPC/SESAM in n/ = 2 
lattice QCD [94]. Lattice results on nucleon GPDs published since then have provided 
important insights into the transverse structure of unpolarized nucleons [95], the lowest 
moments of polarized [96] and tensor GPDs [97], and transverse spin densities of quarks in 
the nucleon [98, 99]. With the exception of several initial studies [100, 101], all previously 
published lattice results on GPDs have been obtained from calculations in a two-flavour 
'heavy pion world' with pion masses in the range of 550 to over 1000 MeV. In Ref. [48] 
previous studies have been improved by presenting a comprehensive analysis of the lowest 
three moments of unpolarized and polarized GPDs in n/ = 2-1-1 lattice QCD with pion 
masses as low as 350 MeV and volumes as large as (3.5 fm)^ (for further reading, see 
also [102, 103, 104, 105, 106, 107, 108, 109]). 

4'2. Double distributions. - The most popular choice when parametrizing GPDs is 
to use a factorized form with a t-dependent part determined by some form factors and 
a t-indcpcndcnt part given in terms of double distributions [71, 72]. Although this as- 
sumption is not strictly valid, it simplifies the QCD evolution considerably because in 
this way the t dependences of quarks and gluons (which mix under evolution) are not 
modified during evolution [110, 111]. The t-independent part H''{x,^) = H'^{x,£,,t = 0) 
is parametrized by a two-component form: 



(4.1) 




where 




1-1/31 



(4.2) 



da6{x- l3-a^) F'^{(3,a). 



-1+1/31 



A similar decomposition is assumed for the spin- flip quark GPD E'^ [16], i.e. 



(4.3) 




where 




1-1/31 



(4.4) 



da5{x- K''{P,a). 



-1+1/31 
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The D-term contribution in Eqs. (4.1) and (4.3) completes the parametrization of 
GPDs, restoring the correct polynomiaUty properties of GPDs [14, 35]. It has a support 
only for \x\ < |^|, so that it is invisible in the forward limit. The _D-tcrm contributes 
to the singlet-quark and gluon distributions and not to the non-singlet distribution. Its 
effect under evolution is at the level of a few percent [112]. 

According to Radyushkin's suggestion [110], the double distributions (DDs) entering 
Eqs. (4.2) and (4.4) can be written as 

(4.5) Fi{p, a) = h{(3, a) q{(3), K<'{p, a) = h{(3, a) e{(3), 

where q{f3) is the forward quark distribution (for the flavour q), while the spin-flip par- 
ton distribution e{(3) cannot be extracted from DIS data and has to be modeled (see, 
e.g. [113]). The profile function h{(3,a) is parametrized as [111] 



22fe-Mr2(6+i) (i_|/3|)2fc-M • 

In Eq. (4.6), the parameter b determines the width of the profile function h{(3,a) and 
characterizes the strength of the ^ dependence of the GPDs. It is a free parameter for 
the valence {hvai) and sea {baea) contributions to GPDs. In such an approach bvai and 
bsea can be used as fit parameters in the extraction of GPDs from hard electroproduc- 
tion observables. The favoured choice is byai — ^sea — 1.0, Corresponding to maximum 
skeweness. With a similar assumption adopted for the gluon distribution one defines 
bgiuon = 2. The limiting case b ^ ao gives ft(/3, a) — > S{a) h{P) and corresponds to the 
^-independent ansatz for the GPD, i.e. H'^{x,^) H'^{xb,^ = 0) = q{xB), as used in 
Refs. [114, 115]. 

Assuming a phenomenological parton distribution one can therefore construct model 
GPDs which include all general constraints and are fiexible enough to allow for a fit of 
the different observables [114, 115, 116, 117, 16, 74]. 

Starting from a similar ansatz in the DGLAP region a parametrization has been 
proposed [118] that satisfies the requirements of polynomiality and positivity conditions 
and is also suitable to study twist-3 effects. Other parametrizations were proposed in 
Ref. [119]. 

Including the D-term the double-distribution representation satisfies polynomiality 
but does not guarantee positivity. An integral representation obeying both polynomiality 
and positivity has been worked out in Ref. [56] on pure mathematical grounds and in 
Ref. [120] from an analysis of simple perturbative graphs for GPDs. Further proposals 
to construct double distributions with dynamical content in connection with the Fock 
expansion of light-cone wave functions can be found in Refs. [121, 122]. 

4'3. The dual parametrization of GPDs. - An alternative way to parametrize GPDs is 
the so-called dual parametrization of GPDs [73] . It is based on the partial wave expansion 
of the GPDs in the t channel [123]. GPDs are then presented as an infinite series of t- 
channel exchanges, which reminds us of the idea of duality in hadron-hadron scattering. 
The applicability of such a procedure is justified by recognizing crossing relations between 
the GPDs of DVCS and the generalized distribution amplitudes of 7*7 hh. 

In fact, the 7*/i jh (DVCS) amplitude at large and small squared momentum 
transfer between the hadrons is related by crossing symmetry to the crossed channel 
amplitude of the 7*7 hh process with a highly virtual photon near threshold, i.e. in 
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the regime where the squared cm. energy is much smaUer than the photon virtuahty 
Q^. At tree level, similar to the DVCS case [8, 9, 10], the 7*7 amplitude can be written as 
a convolution of a hard photon-parton scattering amplitude and a generalized distribution 
amplitude for the soft transition from parton to hadrons [124]. As in the case of DVCS, 
in order to give a nonzero 7*7 hh amplitude the virtual photon must have the same 
transverse hclicity of the real one as a consequence of chiral invariance in the collinear 
hard scattering process. 

The dual parametrization fulfills the polynomiality condition and allows for flexible 
modeling of the t dependence of the GPDs designed primarily for small and medium-size 
values ofa;_B, xb < 0.2 [73]. In addition, crossing is related to the possibility of establish- 
ing dispersion relations between the real and imaginary parts of the elementary DVCS 
amplitude [37, 38] that are helpful to model GPDs. It must be said, however, that math- 
ematically one can not completely restore the complete image of the nucleon contained 
in the GPDs from the knowledge of the elementary DVCS amplitude, because one has to 
deal with an inversion problem. This is a typical problem of tomography which is usually 
solved with the help of Radon transformations [125]. The dual parametrization combined 
with the Radon transformation clarifies which part of the image can be analyzed [39]. 

4"4. Constraints from observables. - The connection between GPDs and Dirac and 

Pauli form factors through the first McUin moment (2.50) is suggesting a simple factorized 
ansatz for the x and t dependences of the GPDs at ^ = 0, as in the VGG model [115, 
114, 116]. Here one assumes 



where Uy{x) and dv{x) are the (unpolarized) distributions for valence up and down 
quarks, respectively. A similar expression is assumed for E"^. The ansatz for H'^ is 



(4.8) H^ix, ^ = 0,t) = Au, 9l{t)lgl{Q), H'^ix, ^ = 0,t) = Ad, 5I WM(0), 



where Auv{x) and Adv{x) are the polarized (helicity) distributions of valence up and 
down quarks, respectively. The quark distributions are taken from parametrizations 
fitting the world DIS data. The ^ dependence is obtained in the VGG model by matching 
this ansatz with the double-distribution approach. 

However, a complete factorization of the x and t dependences of GPDs seems rather 
unrealistic. At small t and small x one can expect Regge behaviour of H'i{x, ^ = 0,t), so 
that one can assume an exponential ansatz: 



(4.9) H'i{x,^ = 0,i) = q,{x) exp[t /"(ar)], E'i{x,i = Q,t) = el{x) exp[tff9(ar)]. 



(4.7) 



H''{x,^ = Q,t)=u,{x)F^{t)/2, H\x,^ = Q,t)=d,{x)Ft{t) 



where e1{x) is just the forward limit of E'i{x,^ = 0,t) normalized as 



(4.10) 




and the fimctions f'^{x) and g'^{x) parametrize how the profile of the quark distribution 
in the impact parameter plane changes with x. Specific forms have been proposed for 
them in Refs. [74, 75] in order to be consistent with dominance of the leading meson 
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trajectories known from Regge phenomenology and to satisfy the Drell-Yan-West rela- 
tion [126] between the large-x power behaviour of parton distributions at t = and the 
largc-f power behaviour of the associated clastic form factors. With only four parameters 
in Ref. [74] and with much more flexibility in Ref. [75] a rather good fit to elastic nucleon 
electromagnetic form factors is obtained. Considerable ambiguities have been found in 
determining e1{x) since the forward limit of E'^(x,^,t) is not known. The form factor 
data disfavour an identical shape in x and t of E^{x,S^ = 0,t) and E'^{x,£, = 0,t) [75]. 

A physically motivated parametrization for the unpolarized quark GPDs H and E 
has been recently proposed [76] using a combination of constraints from simultaneous 
fits of the experimental data on both the nucleon elastic form factors, the DIS structure 
functions in the non-singlet sector and the lattice calculations of Mellin moments with 
n> 1. The method is thus able to produce ^- and f-dependent — H'^ and E'^ — E"^ 
satisfying the polynomiality condition by construction. 

4'5. The chiral quark- soliton model. - The chiral quark-soliton model (xQSM) [127] 
is based on the principles of chiral symmetry breaking and the limit of a large number of 
colours Nc- The underlying effective relativistic quantum field theory was derived from 
the instanton model of the QCD vacuum [128], which provides a mechanism of dynamical 
chiral symmetry breaking, and is valid at low energies below a scale of about p~J^ « 600 
MeV, where pav is the average instanton size. Thus the xQSM can be considered as a 
realization of the idea that in the large- TV^ limit the nucleon can be viewed as a classical 
soliton of the pion field. 

In the leading order of the large- A'^c limit the pion field is static, and one can determine 
the spectrum of the effective one-particle Hamiltonian of the theory: 

(4.11) Hes\n) = En\n), H,s = -il°l''dk + 7°MC/^^ 

where W^'' — exp(i75T°7r") with U = exp(iT"7r°) denoting the SU{2) chiral pion field. 
The spectrum consists of an upper and a lower Dirac continuum, which are distorted by 
the pion field as compared to continua of the free Dirac-Hamiltonian, 

(4.12) Ho\no) = -E„Jno), Ho = -ztV^/c +7°M, 

and of a discrete bound state level of energy E\cv By occupying the discrete level and 
the states of lower continuum each by Nc quarks in an anti-symmetric colour state, one 
obtains a state with unity baryon number. The soliton energy Egoi is a functional of the 
pion field. 



(4.13) E,oi[U]=Nc 



+ ^ (En — Ena) 



Minimization of i?soi[t^] determines the self consistent solitonic pion field Uc- The nucleon 
mass Mjv is given by Eso\ [Uc] ■ Quantum numbers of the baryon - like momentum, spin 
and isospin - are described by considering zero modes of the soliton. Corrections inl/Nc 
can be included by considering time-dependent pion field configurations. The results of 
the xQSM respect all general counting rules of the large- A^c phenomenology. 

In the xQSM one can evaluate in a parameter-free way nucleon matrix elements of 
QCD quark bilinear operators as 



(4.14) {N',p'\%{zi)ri,g{z2)\N,p) 
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f dPX e*(p'-p)-^ - X)r$„(z2 - X) e' 



where for zi ^ Z2 the insertion of the gauge Knk is understood on the left hand side. The 
dots in Eq. (4.14) denote terms subleading in the l/N^ expansion. In Eq. (4.14) F is some 
Dirac and flavour matrix, A^vw' constant depending on F, the spin and flavour quantum 
numbers of the nucleon state | A^) = jS'a, T3), and <I>„(a;) = {x\n) are the coordinate-space 
wave-functions of the single quark states \n) defined in Eq. (4.11). The sum in Eq. (4.14) 
goes over occupied levels n (i.e. n with En < i^iev), and vacuum subtraction is implied 
for En < E\ev as in Eq. (4.13). 

Many static nuclconic obscrvablcs, like magnetic moments, electric polarizabilitics, 
axial properties, have been computed in the xQSM in the way sketched in Eq. (4.14). 
The results were found in good agreement with data (see Ref. [129] for a review). In 
particular the model describes data on electromagnetic form factors up to |t| ^ 1 GeV^ 
within (10-30)% [130, 129]. In Ref. [131] it has been demonstrated that the model can 
be applied to the description of twist-2 quark and anti-quark distribution functions of 
the nucleon. The consistency of the approach has been shown by giving proofs that 
the model expressions satisfy all general requirements of QCD [131]. The distribution 
functions computed in the xQSM [131, 132, 133] (and also [134, 135, 136]) refer to a 
low normalization scale of around 600 MeV, and agree with available parametrizations 
performed at comparably low scales [137] within (10-30)%. 

The xQSM lacks explicit gluon degrees of freedom. Therefore one deals with quark 
GPDs only. Different flavour combinations of the GPDs exhibit different behaviour in 
the large- iVc limit [16] 

(4.15a) (F" + H''){x, e, t) = N^fiu, v, i), (£;" - E''){x, t) = Nlf{u, v, t), 
(4.15&) {H^-H''){x,i,t)=N,f{u,v,t), iE^+E^)ix,^,t)=N^f{u,v,t), 
(4.15c) {H"-H''){x,i,t)=N^f{u,v,t), {E^ - E^){x, ^,t) = fiu,v,t), 
(4.15rf) {H^ + H^){x,^,t)^N,fiu,v,t), iE"+E^){x,^,t)=Nlf{u,v,t). 



The functions f{u,v,t) are stable in the large- A^c limit for fixed values of the 0{N^) 
variables u = NcX, v = Nc^ and t, and of course different for the different GPDs. 

The model expressions for the leading GPDs in Eqs. (4.15a) and (4.15c) were derived 
in Refs. [81, 82], respectively, where it was also demonstrated that they are correctly nor- 
malized to the corresponding electroweak form factors. In Ref. [138] the expressions for 
the leading (chiral-even) GPDs have been proved to also satisfy the polynomiality prop)- 
erty. Numerical estimates show the relevant role of the Dirac continuum in determining 
the shape of the GPDs (Fig. 3) . The forward limit of the unpolarized spin-flip isoscalar 
and isovcctor GPDs, [E"^ + E'^){x,^,t) and (£^" - E'^){x,^,t), have been computed in 
Refs. [83] and [84], respectively. 

4'6. Constituent- quark models. - The modeling of GPDs in constituent quark models 
(CQMs) relies on the basic assumption that there exists a scale Qq where the short 
range (perturbative) part of the interaction is negligible and, therefore, the glue and 
sea are suppressed, while a long range (confining) part of the interaction produces a 
proton composed by (three) valence quarks, mainly [139]. Jaffe and Ross [140] proposed 
to ascribe the quark model calculations of matrix elements to that hadronic scale Qq, 
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Fig. 3. - The isosinglet distribution H{x,i,A'^) = H"{x,i,A'^) +H''{x,^,A^) in the forward 
hmit, = and ^ = 0, (left panel) and for A5_ = -A^ - ^^M^ = and ^ = 0.3 (right 
panel). Dashed lines for the contribution from the discrete level; dashed-dotted lines for the 
contribution from the Dirac continuum; solid lines for the total contribution [81]. 



typically Qq = 0.2 — 0.4 GeV^. In this way, quark models summarizing a great deal 
of liadronic properties may substitute for low-energy parametrizations, while evolution 
to larger momentum is dictated by perturbative QCD [90]. CQMs were applied 
in the calculation of unpolarized and polarized GPDs in Ref. [87] in terms of overlap of 
Schrodinger wave functions for three constituent quarks obtained within the instant- form 
quantization and for nonrelativistic kinematics. A fully covariant formalism has been 
adopted in Refs. [88, 89, 91] for the calculation of the GPDs both in chiral-even and chiral- 
odd sectors within the framework of light-cone quantization, using the representation of 
GPDs in terms of overlaps of light-cone wave functions (LCWFs) describing the A^-parton 
composition of a hadronic state [57, 141]. The valence-quark contribution is obtained 
by specializing to the case N = 3. The corresponding valence-quark component of the 
LCWF can be related to a wave function obtained in the (canonical) instant form through 
appropriate Melosh rotations. Therefore, one can link LCWFs to wave functions derived 
in CQMs where they are obtained as eigenfunctions of the nucleon Hamiltonian in the 
instant-form dynamics. Of course, this link is useful in the kinematic range where valence 
quark degrees of freedom are effective. In this region GPDs exhibit the exact forward 
limit, reproducing the parton distribution with the correct support and automatically 
fulfilling the particle number and momentum sum rules. In addition, from the LCWF 
overlap representation of the PDFs for the valence-quark contribution, it is possible to 
derive the following relations valid at the hadronic scale Qq [91]: 

(4.16) 2Atu{x) = Au{x) + lu{x), 2ATd{x) = Ad{x) - ^d{x). 

These relations are compatible with the Soffer inequality [142] and generalize to the case 
of parton distributions the results obtained in Refs. [143, 144] for the axial (Aq) and 
tensor {6q) charges: 



(4.17) 



25q = Ag + Ag-NR, 
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where A^nr is the axial charge in the nonrelativistic hmit, i.e. 

(4.18) A^NR = (|^.,l/2 - i^r,-l/2) • 

Using LCWFs derived from the relativistic CQM of Ref. [145] the resulting tensor 
charges are 5u = 1.16 and 6d = —0.29 [146], whereas those predicted by QCD sum 
rules [147] are Su = 1.33 ± 0.53 and Sd = 0.04 ± —0.02. When evolved in leading-order 
QCD from the intrinsic scale of the model to — 10 GeV^ they become du = 0.79 and 
Sd = —0.20, within the range of values calculated in the different models considered in 
Ref. [148] and in fair agreement with lattice QCD calculations [106]. 

With the same model the sum rule (2.69) for the angular momentum carried by trans- 
versely polarized quarks gives quite substantial contributions for up and down quarks, 
i.e. [146] 

(4.19) {6''J^) = 0.54, (5^ J|) = 0.37, 

and for the quantity defined in Eq. (2.53) one finds Kj. = 3.98 and Kj. = 2.60. The 
same positive sign of /t^ for both up and down quarks would imply [40] a negative Boer- 
Mulders function h^'^ describing the asymmetry of the transverse momentum of quarks 
perpendicular to the quark spin in an unpolarized nucleon [149]. 

4'7. The meson-cloud model. - The pion cloud of the nucleon associated with chiral- 
symmetry breaking was first discussed in the DIS context by Feynman [150] and Sul- 
livan [151]. As realized by Thomas [152], it can give an explanation of the flavour- 
symmetry violation in the sea-quark distributions of the nucleon, thus naturally ac- 
counting for the excess of d (anti)quarks over u (anti)quarks as observed through the 
violation of the Gottfried sum rule [153, 154, 155, 156]. The connection between the 
meson-cloud model and the chiral properties of GPDs was also studied by investigating 
the nonanalytic behaviour of the d — u distribution in [157, 158]. Although the nu- 
cleon's nonpcrturbativc antiquark sea cannot be ascribed entirely to its virtual meson 
cloud [159], the meson-cloud model for the physical nucleon has a long and successful 
history in explaining properties such as form factors [160, 161, 162] and parton distribu- 
tions [163, 164, 165]. It has been revisited and applied for the first time to study quark 
GPDs in [92]. The meson cloud was also taken into account to study the unpolarized 
gluon distribution in Ref. [166]. 

The basic assumption of the meson-cloud model is that the physical nucleon N is made 
of a bare nucleon N dressed by the surrounding meson cloud, so that the state of the 
physical nucleon can be decomposed according to the meson-baryon Fock-state expansion 
as a superposition of a barc-nucleon state and states containing virtual mesons associated 
with recoiling baryons. This state, with four-momentum = {pJf,P^,PN±) = {Pn^Pn) 
and helicity Ajv, is an eigenstate of the light-cone Hamiltonian 

(4.20) Hlc = Yl [^oiQ) + H^iQ) + Hi{N, BM)] . 

B,M 

In Eq. (4.20) -ff(f (<?) stands for the effective-QCD Hamiltonian which governs the consti- 
tuent-quark dynamics, and leads to the confinement of three quarks in a baryon state; 
analogously, HQ^{q) describes the quark interaction in a meson state, Hi{N, BM) is the 
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nucleon-baryon-meson interaction, and the sum is over all the possible baryon and meson 
configurations in which the nucleon can virtually fluctuate. 

In a pcrturbativc treatment of the meson effects, the Fock-space expansion is trun- 
cated to the one-meson components, and the nucleon wave function is expanded in terms 
of eigenstates of the bare Hamiltonian Hq = HQ{q) + HQ^{q). The corresponding state 
of the physical nucleon \N) can be written as 

(4.21) \pN,\N\N)=\fZ\pNAN\N) 

^ r AyA^k^ 1 ^ A„ (N,BM) ^ ■ 

X \yp+, k± + ypi^^, A'; B)\{1- y)p+, -k± + {1 - y)pM±, A"; M), 

where the function (j)^7y^'''^^^\y,kj_) defines the probability amplitude for a nucleon 
with helicity Aat to fluctuate into a virtual BM system with the baryon having helicity 
A', longitudinal momentum fraction y and transverse momentum fej_, and the meson 
having helicity A", longitudinal momentum fraction 1 — y and transverse momentum 
—k±. The normalization constant Z in Eq. (4.21) gives the probability that the nucleon 
is a bare core state and ensures the correct normalization of the nucleon wave function: 

(4.22) {p'+,p'^,\';H\p+,p^X-H) = 2{27rfp+S{p'+ - p+)S^^\p'^ - p^)Sxx' ■ 

In the region ^ < a; < 1, where a quark is taken out of the nucleon with momentum 
fraction .t + ^ and reinserted back with momentum fraction x — ^, the virtual photon 
can hit either the bare nucleon N or one of the higher Fock states. As a consequence, 
the matrix elements defining the GPDs in Eq. (2.22a) can be written as the sum of two 
contributions, 

where ^^-^are contribution from the bare proton, described in terms of Fock states 

with three valence quarks. This term can be calculated in the light-front overlap represen- 
tation derived in Ref. [57], and applied to the case of three valence quarks in Refs. [88, 89], 
where one can also find explicit expressions in terms of bare-nucleon LCWFs derived in 
a CQM. The 5A'i term in Eq. (4.23) is the contribution from the BM Fock components 
of the nucleon state, corresponding to five-parton configurations. This term can further 
be split into two contributions, with the active quark belonging either to the baryon 
(M«/(^'^)^) or with the active quark inside the meson (M«/(^'^)^), i.e. 

(4.24) 5AI, = y V 5^«/(^'f)^. 

B,B',M M',M,B 

The term is schematically represented in Fig. 4a. It is explicitly given by 

1 A.-i^ r A2i 



(4-25)<^^.^,,.„,(x,^,A^)-^2L.J^ J 2(27r)3^AV,A.(y-^'y^'^^) 
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Fig. 4. - Deeply virtual Compton scattering from the virtual (a) baryon and (b) meson compo- 
nents of a dressed nucleon. 

where ks = {ks + kB>)/'2. is the average momentum of the baryon in the initial and final 
state, and ys = "k^/P^ is the average fraction of the plus-momentum of the nucleon car- 
ried by the baryon. The term ^v^i^a^^ -'^l- (4-25) is the matrix element of Eq. (2.22a) 
evaluated between the initial {B) and final {B') baryon state of the baryon- meson com- 
ponent of the nucleon and can explicitly be calculated in terms of overlap of bare-baryon 
LCWFs. 

Analogously, the contribution from the meson in the BM fluctuation (Fig. 4b) is 
calculated by interchanging the role of the baryon and meson substates. It reads 

(4-26)M,,^,,,„^(x,e,A^)-^2.J^ 2(27r).^^^V,A. (-^,-^,A^) 

with ^v^^^A^T^ matrix element between the initial and final meson participating 
to the interaction process. As a result, the contribution from the BM components in 
Eqs. (4.25) and (4.26) are obtained by folding the matrix elements of the baryon and 
meson constituents with the probability amplitudes describing the distributions of these 
constituents in the dressed initial and final nucleon. We also note that the sum in 
Eq. (4.24) over all the possible BM configurations includes contributions from both the 
diagonal matrix elements with the same hadrons in the initial and final state {B' = B 
and M' = M in Eq. (4.25) and (4.26), respectively), and the matrix elements involving 
the transition between different hadron states (i.e. the terms with B' ^ B and M' ^ M 
in Eq. (4.25) and (4.26), respectively). 

In the region —l<x< where emission and reabsorption refer to an antiquark 
with momentum fraction —{x + £,) and —(a; — ^), respectively, the only nonvanishing 
contribution comes from the active antiquark in the meson substate of the BM Fock 
component of the nucleon wave function, i.e. 

(4-27) ^I,,„(x,^,A,)= ^ Mf/r:-)-(.,^,A,), 

M,M',B 

where (JAy ^^^'^^''^ is given by the convolution formula (4.26) taken for negative value of 

X. 
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In the region — ^ < a; < ^ the scattering amphtude describes the emission of a quark- 
antiquark pair from the initial nucleon, so that in the Fock-state decomposition the parton 
content of the initial nucleon must differ from the final one by such a qq pair. In the 
meson-cloud model, the final state is described by the three-valence quark configuration 
of the bare nucleon and the initial state is given by the five-parton component of the 
BM substatc, with the active quark belonging to the baryon and the active antiquark 
to the meson. The contribution with both the active quark and antiquark belonging to 
the meson substate of the initial nucleon is vanishing by orthogonality between the final 
bare- nucleon state and the initial baryon state [92] . 

This convolution formalism is applied here to calculate the polarized and unpolarized 
GPDs, taking into account the meson-cloud contribution corresponding to n, p, and w, 
with the accompanying baryon in the BM component of the dressed nucleon being a 
nucleon or a A. For the bare-hadron states we use the LCWFs specified in Ref. [162], 
which are able to give an overall good description of the nucleon electroweak form factors. 

In Fig. 5 the spin-averaged H'^ and the helicity-flip E'^ GPDs arc plotted as function 
of X at different fixed values of ^ and t together with the separated contributions from 
the bare proton (dashed-dotted lines), the baryon (dashed lines) and the meson (dotted 
lines) in the baryon-mcson fluctuation. All these contributions add up incoherently to 
give the total result (full curves). The bare proton is always positive within its support 
(0 < X < 1) with the exception of E"^ for which it is negative. The same behaviour 
characterizes the baryon contribution from the baryon-meson fluctuation that is also 
limited to the range < a; < 1, consistently with the assumption that the only active 
degrees of freedom for such a baryon are the valence quarks. Both contributions vanish at 
the end points of their support. The contribution from the mesons in the higher Fock state 
components extend all over the full range — 1 < x < 1, and it is determined at positive 
(negative) x by the quark (antiquark) residing in the mesons. The resulting effect of the 
meson cloud is thus to add a contribution for negative x and to increase the magnitude 
of the GPDs for positive x with respect to the case of the bare proton. In particular, 
for positive and small x the meson-cloud contribution as a whole is comparable to that 
of the bare proton, confirming the important role of the sea at small x found within 
the chiral quark-soliton model [81, 83, 84]. Although we adopted instant-form hadron 
wave function which are SU(6) symmetric, the hclicity flip GPDs arc different from 
zero because of the orbital angular momentum generated by the Melosh rotations in the 
transformation from the canonical spin to the light-front spin. Therefore the faster fall off 
of -E"? with respect to ff"? for .t ^ 1 is a consequence of the decreasing role of the Melosh 
transform to generate angular momentum in E'^ with increasing quark momentum. From 
the first moment of E'^ at t = and ^ = we obtain the following values for the quark 
anomalous magnetic moment: k" = 1.94 and k"^ = —1.73, corresponding to = 1.87 
and At" = —1.80. These values include a contribution from the higher Fock states of the 
order of 15%, which is essential to bring the values of pretty close to the experimental 
values kP = 1.793 and k" = -1.913. 

Although we use a covariant formalism based on the LCWF overlap representation of 
GPDs, one notices from Fig. 5 that the first moment of H'^ and E"^ is not ^-independent 
as required by Lorentz invariance. This is a consequence of the model calculation where 
we introduced phenomenological LCWFs which violate Lorentz symmetry. As outlined in 
Ref. [167] , the functional form of LGWFs is dictated by the underlying Lorentz symmetry, 
i.e. the longitudinal and transverse variables are tied to each other in a certain but non 
apparent manner. Hence, phenomenological LCWFs may result in a violation of the 
GPDs polynomiality property. 
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Fig. 5. - The spin averaged {H'', three upper panels) and the helicity flip {E'', three lower 
panels) generalized parton distributions for the u (left panels) and d (right panels) flavours, at 
fixed values of ^ and t as indicated. Dashed-dotted curves: bare contribution. Dashed curves: 
contribution from the active baryon. Dotted curves: contribution from the active meson. Solid 
curves: total result. 



With a nonvanishing ^ one can explore the ERBL region with |a;| < ^ where in our 
model only transition amplitudes between the bare-proton and baryon-meson components 




Fig. 6. - The helicity dependent spin averaged {H'', three upper panels) and the hehcity flip {E'', 
three lower panels) generalized parton distributions for the u (left panels) and d (right panels) 

flavours, at fixed values of ^ and t as indicated. Dashcd-dotted curves: bare contribution. 
Dashed curves: contribution from the active baryon. Dotted curves: contribution from the 
active meson. Solid curves: total result. 



are contributing. We see that the GPDs in the ERBL region are rather regular functions 
over the whole range, with zeros at the endpoints x = This result is quite different 
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from the oscillatory behaviour predicted by the chiral quark-soliton model [81] where the 
valence contribution of the discrete level is a smooth function extending into the ERBL 
region and adding to the sea contribution. Here this is forbidden because the support of 
the valence contribution is limited to the DGLAP region. 

The results for the bare-nucleon and BM-component contributions to the helicity- 
dcpcndcnt GPDs are shown in Fig. 6. Comparing the results at ^ = and different t 
values, we see that both and H'^ exhibit a small t-dependence. At constant t = —0.2 
GeV^, the spin-averaged GPDs H have a rather weak ^-dependence, at variance with 
the behaviour of the spin-flip GPDs E. Due to the opposite sign of i7" and H'^ in all 
kinematic conditions, their difference is positive and peaked at a value of x comparable 
to the results obtained in the chiral quark-soliton model in the leading order of the l/N^ 
expansion [82]. Also i?" and E'^ have opposite sign as functions of .t, and at small a; > 
exhibit a negative and positive peak, respectively, with a sign change at intermediate 
value of X. The contribution in the ERBL region from the interference between BM 
highcr-Fock states and the barc-nuclcon state is more pronounced in the case of the 
spin-flip GPDs. The isovector combination E"^ — E'^ in the ERBL region has the same 
negative sign as in quark-soliton model, but it is much smaller in absolute value [82]. We 
also note that in our calculation the pion-pole contribution, giving a large contribution 
in the ERBL region, is not taken into account. 

At negative x, the only contribution is from the antiquark of the meson substates. 
In the case of the hclicity-dependent GPDs, the non- vanishing contributions from the 
diagonal terms with the same meson in the initial and final states are from the p and oj 
mesons, while the pion can enter only in the interference terms with the w and p. Since the 
probability to find the p or the oj meson in the dressed nucleon is much smaller than that 
for the pion [162], we find that the meson contribution to the helicity-dependent GPDs 
is much smaller than in the case of the unpolarized GPDs, and it is almost vanishing at 
all values of x. 



Table II. - Contributions to the angular momentum sum rule (2.62) from the bare nucleon 
(second column), the meson cloud at positive x (third column), the meson cloud at negative x 
(fourth column), and the total result (last column). 



bare nucleon 



BM {x > 0) 



BM {x < 0) 



TOT 



^2,o(0) 

^lo(O) 

B?.o(0) 
BioiO) 



0.61 
0.30 
0.43 
-0.43 
0.90 
-0.23 
0.07 
0.05 



0.047 
0.032 
0.027 
-0.040 
0.062 
-0.014 
0.006 
0.003 



0.0055 
0.0086 
-0.0026 
0.015 
-0.00087 
0.00018 
0.002 
0.012 



0.66 
0.34 
0.45 
-0.45 
0.96 
-0.24 
0.078 
0.065 
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The different GPD moments entering the angular momentum sum rule (2.62) are 
presented in Table II, with the separate contributions from the bare nucleon and the 
higher Fock-state components at a; > and x < 0. In all cases the contribution from 
the BM components of the dressed nucleon are within a few percent of the dominant 
bare- nucleon contribution. We note that S" is much larger than E**, while the angular 
momentum contribution is small and of the same amount for both flavour quarks, despite 
the fact that there are less down quarks in the proton. Consistently with the fact that 
in our model we do not have gluons and therefore the nucleon momentum is saturated 
by quarks and antiquarks, we flnd _B| q(0) = 0. 

These results are all at the hadronic scale of the model, and evolution should be 
applied before comparing with other model calculations and experiments referring at 
higher scales. 

5. — Applications to observables 

5T. Electroweak form factors. - Electroweak form factors of the nucleon have played 
a privileged role in the investigation of the nucleon structure and a variety of models 
have been devised to accoimt for data (for recent reviews, sec [168, 169, 170, 171, 172]). 
In the present context they are considered as first Mellin moments of GPDs. As such 
they have been used as constraints in building phenomenological GPDs as discussed in 
sect. 4.4. Conversely, model GPDs can be tested in a comparison with the observed form 
factors as done in a recent investigation based on the meson-cloud model [162]. 

According to the analysis of Ref. [173] a pronounced bump structure in the neutron 
electric charge form factor G'^(Q^) (and a dip in the other nucleon form factors) aroimd 

= 0.2 — 0.3 GeV^ can be appreciated and interpreted as a signature of a very long- 
range contribution of the pion cloud surrounding the bare-nucleon core and extending 
out to 2 fm (^). This behaviour has been reproduced in the Lorentz covariant quark 
model of Ref. [174]. It must be said, however, that from dispersion relation analysis [175] 
the pion cloud should peak much more inside the nucleon, at ~ 0.3 fm, and the desired 
bump-dip structure of Ref. [173] can only be achieved at the cost of low- mass poles close 
to the u) mass in the isoscalar channel and to the three-pion threshold in the isovector 
channel [176]. 

As a matter of fact, while giving an overall good description of all the electroweak 
form factors, the meson-cloud model is unable to reproduce such a bump/dip structure 
(Fig. 7) . In any case the contribution of the meson cloud is smooth and only significant 
below = 0.5 GeV^. For the neutron electric form factor it is essential to also include 
a small percentage (1 — 2%) of mixed-symmetry 5'-wave momentum component in the 
three-quark wave function, in agreement with earlier findings [177] and only slightly 
affecting the other form factors. 

Charge and magnetization densities are encoded in the electromagnetic form factors, 
but their radial distribution cannot be extracted without ambiguities [179]. Neglecting 
relativistic corrections, in the Breit frame they are given by the Fourier-Bessel transform 
of the nucleon electromagnetic Sachs form factors. In the meson-cloud model the meson 
cloud manifests itself as a slight extension of the radial distribution up to ~ 2 fm. In 
the neutron case the resulting charge distribution shows a positive core surrounded by a 
negative surface charge pushed outwards by the meson cloud and peaking at ^ 0.8 fm, 



(^) Throughout this section the notation = 



—t is used for the four-momentum transfer. 
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Fig. 7. — The four nucleon electromagnetic form factors compared with the world data considered 
in the analysis of Ref. [173] and the recent JLab data [178] using G% = {fjFG%/Gl^)/(l + 
/Q.llGeV^Y (open squares). Long-dashed curves for the contribution of the meson cloud; 
dashed curves for the valcnco-quark contribution with SU(6) instant-form nucleon wave function; 
solid curves for the sum of the two contributions; dashed-dotted curves for the total result with 
1% mixed-symmetry S"-state in the bare nucleon wave function (taken from Ref. [162]). 



in agreement with the analysis of Ref. [179] and the expectation based on the picture of 
a hadron's periphery caused by the pion cloud [180]. 

In contrast, the density q{b) of partons of type q in the transverse (impact parameter) 
plane with respect to the direction of the three- momentum transfer can be unambiguously 
determined by taking the two-dimensional Fourier transform of the Dirac form factor F^: 

1 /"^ 

(5.1) q{h) = —j^ dQQJo{Qb)F?iQ^), 

where Jo is a cylindrical Bessel function [181]. Assuming that only up and down quarks 
are in the nucleon and invoking isospin symmetry, from the relations 

(5.2) pPi^t>) = ^u{b)-^d{b), p»(6) = -^«(6) + ^d(6) 

one derives the proton and neutron transverse charge distributions. In Eq. (5.2), u{b) 

and d{b) arc the up and down quarks distributions in the proton, which arc the same as 
for down and up quarks in the neutron, respectively (with normalization / d^bu{b) = 
J d^bdib) = 1). The corresponding charge densities for the proton and the neutron are 
plotted in the upper panels of Fig. 8. As in the phenomenological analyses of Refs. [181, 
182] the densities are concentrated at low values of b with a positive peak for the proton 
and a negative peak for the neutron. While the negative tail of the neutron distribution 
has the well known interpretation in terms of the pion cloud, the negative charge density 
near the origin appears to be mysterious. An intuitive understanding of this result has 
been recently given in Ref. [183] , suggesting that the up quark in the neutron has a larger 
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0.5 1 1.5 2 0.5 1 1.5 2 

b (fm) b (fm) 

Fig. 8. - Upper panels: The proton (left) and neutron (right) charge densities as a function of 
the impact parameter b. Lower panels: The up (left) and down (right) quark transverse charge 
densities in the proton. Solid curves for a SU(6)-symmctric instant-form wave function, dashed 
curves with mixed-symmetry components (taken from Ref. [162]). 




0.5 1 1.5 2 0.5 1 1.5 2 

b (fm) b (fm) 



Fig. 9. - Transverse distribution of up and down quarks in a longitudinally polarized proton 
as a function of the impact parameter b. Left panels: the axial contributions Au and Ad for 
up and down quarks, respectively. Right panels: total contribution for quarks polarized in the 
longitudinal direction, either parallel (solid lines) or antiparallel (dashed lines) to the proton 
helicity (taken from [162]). 



p-wave component than down quarks, being therefore suppressed at the origin as shown 
in the lower panels of Fig. 8. 
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By Fourier transforming the quark contributions to the Dirac and axial form factors 
one has access to the probabiUty A, Aat) to find a quark with transverse position b 
and light-cone helicity A (= ±1) in the nucleon with longitudinal polarization Ajv (= ±1), 
i.e. 

(5.3) p^ib, A, \N) = \j [F^iQ^ = qi) + X\n9\{Q^ = qD] e'"-'" 

= 1 [q^b] + X\NAq{b)] , 

where q{b) was already defined in Eq. (5.1) and Aq{b) is the Fourier transform of g^((5^), 
normalized as / cPbAu{b) = 0.96 and / (fbAd{b) = —0.24 (see Table II). Assuming a 
positive proton helicity (Ajv = 1) the resulting probability is shown in the right panels of 
Fig. 9. The axial contributions Au{b) and Ad{b) for up and down quarks (left panels), 
respectively, have opposite sign. When suitably combined with the corresponding trans- 
verse distributions u{b) and d{b) the positive helicity up quarks in the proton are found 
to be preferentially aligned with the proton helicity, while the opposite occurs for down 
quarks. 

5'2. Form factors of the energy-momentum tensor. - The nucleon form factors of the 
energy-momentum tensor were subject to modest interest in the literature for a long time. 
After some pioneering analysis of the mass structure of the nucleon in terms of quark and 
gluon contributions to matrix elements of the energy- momentum tensor [184], only very 
recently, in connection with the advent of GPDs accessible in hard exclusive processes, 
their importance has been appreciated. Generalized form factors of the quark part of the 
energy-momentum tensor have been calculated in lattice simulations in Refs [185, 186] 
and more recently by the LHPC [94, 95, 101, 48, 105] and QCDSF [93, 97, 187, 188] 
collaborations. Quite recently they have also been studied in the xQSM [44, 189, 190] 
and in the Skyrme model [191]. 

In lattice simulations a comprehensive study of the lowest moments of the generalized 
form factors is possible at present for pion masses as low as 350 MeV and volumes as large 
as (3.5 fm)"^ [48], providing results at the scale /U^ = 4 GeV^ in the MS renormalization 
scheme. 

The t dependence of the different Mellin moments can be fitted to a dipolc form, 
with a clear fiattcning of the t slope with increasing order of moment [95, 103]. The 
absolute values in the isovector and isosinglet channels are in qualitative agreement with 
the predictions from large Nc counting rules for the unpolarized generalized form factors, 
i.e. 

(5.4) 1^2"+'^] ~ » lC»o'| ~ iVe, 

whereas in the polarized case counting rules are only partially satisfied [48] . 

Extrapolating to the physical pion mass requires a combination of full QCD lattice 
calculations and chiral perturbation theory (ChPT). Significant progress has been made 
in several approaches to ChPT like heavy-baryon ChPT (HBChPT) [192, 193, 194], 
covariant ChPT in the baryon sector (BChPT) [195], self-consistently improved one- 
loop BChPT [196], and BChPT with finite regulators [197, 198]. As an example, in 
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0.2 0.4 0.6 O.S 0.2 0.4 C.6 0.8 

m^lGeVl im,[GeV| 

Fig. 10. - Fit of the isovector moments A^^'^iO) = {x)u-d and 0*^(0) = (Ax>„_d in the 0{p^) 
covariant BChPT [195] to the LHPC lattice data of Ref. [48] (solid line). The bands shown 
indicate estimate of higher order possible corrections. The lattice results refer to a scale /i'^ = 4 
GeV^ (taken from [195]). 



Fig. 10 the isovector moments g''(0) and q''(0) are shown as a function of the pion 
mass obtained in lattice calculations [48] and in the covariant BChPT at O^p^) [195]. 
The chiral curvature in both observables naturally bends down to the phenomenological 
value for lighter quark masses, leading to a very satisfactory extrapolation curve. 

The generalized form factors ^2,0 ('^ — ^) ^"^"^ ^2o(* = ^) enable us to compute the 
total angular momentum contribution to the nucleon spin J'^ according to Ji's sum rule 
(2.62). In turn, the quark spin is given by ^AE? = iif_o(* = 0) = ^ / H«(a;, 0, 0), so 
that the orbital angular momentum can be derived from lattice simulations versus the 
pion mass m,r- Two remarkable features are found. The first is that the magnitude of 
the orbital angular momentum contributions of the up and down quarks are separately 
sizable, i" « —U^ « 0.30, yet they cancel nearly completely at all pion masses, m 
[48, 199]. The second is the close cancellation between the orbital and spin contributions 
of the down quarks for all pion masses, J"* « [48, 199]. Therefore these results indicate 
that the total angular momentum of quarks in the nucleon is of the same size as the 
quark spin contribution, while the orbital angular momentum is consistent with zero. 

In the Breit frame characterized by A° = one defines the static energy-momentum 
tensor for quarks (and analogously for gluons) [44], 

(5.5) e^''(r, ^) = -L exp(iA • r) (p', 5'|0,^''(O)|p, S), 

with the initial and final polarization vectors of the nucleon S and S' defined such that 
they are equal to (0, s) in the respective rest-frame, where the unit vector s denotes the 
quantization axis for the nucleon spin. 

The components of 6q,ofe('",'S) and e''^'^rjQqfik{r,s) correspond respectively to the 
distribution of quark momentum and quark angular momentum inside the nucleon. The 
components of {Qq.ik — ^Sik'^q.u)ii'i characterize the spatial distribution of 'shear 
forces' experienced by quarks inside the nucleon. The respective form factors are related 
to e^=^(r,s) by 

(5.6a) Jq{t) + ^Jq'{t)= /dVe-''-^£'^"'=s,r,e,,ofe(r,s), 
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(5.66) dUt) + I dfit) + ^ dfit) = / dVe— ^ e,,,(r) (rV^' - ^ 6'^^ , 

(5.6c) M2{t) - ^ (^M2{t) - 2J{t) + ^di(i)) = ]\^/ rfVe-''-^eoo(r,s), 

where the prime denotes derivative with respect to t. For a spin-i particle only the 
0°^-components are sensitive to the polarization vector. Eq. (5.6c) holds for the sum 
Goo = E,S,,oo + Bg,oo with M2{t) = J2qMl(t) + Ml{t) and J(i) and di(t) defined 
analogously, but not for the separate quark and gluon contributions, since otherwise the 
form factor c{t) would not cancel out. 

The form factors M|'^(t), Jq,g{t) and df^{t) depend on the renormalization scale /i. 
Their quark+gluon sums, however, are scale independent form factors. For m,r 7^ 
within the framework of the xQSM, M2{t), J(t) and di(t) can well be approximated by 
dipole fits [44]. J{t) exhibits a similar t dependence as the electric Sachs form factor 
G£;(f). However, M2{t) falls off with increasing \t\ slower than GE{t), while di{t) shows 
a faster fall off. These results indicate that factorizing the t dependence of GPDs is quite 
a rough approximation in the xQSM. 

At t = M2{t), J{t) and di{t) satisfy the constraints, 

(5.7a) M2(0) = ^ Jd^'r Qooir, s) = 1, 

(5.76) J(0) = J d^r e'^^ Si rj Ookir, s) = ^ 



2' 



(5.7c) di(0) = ^ / d^r Oi.ir) ( rV^' - — 6'^ ] = d^. 



Eqs. (5.7a) and (5.76) simply mean that in the rest frame the total energy of the nucleon 

is equal to its mass and the spin of the nucleon is ^. The value of di is not known a priori 
and must be determined experimentally. Remarkably, di determines the behaviour of 
the D-term (and thus the unpolarized GPDs) in the asymptotic limit of renormalization 
scale fi ^ oo [16]. 

The form factor di{t) is connected with the distribution of pressure and shear forces 
experienced by the partons in the nucleon [200], as can be verified by recalling that Qij{r) 
is the static stress tensor which (for spin and ^ particles) can be decomposed as 



(5.8) eij{r)=s{r) Af--dij)+p{r)6., 



nrj_ _ 1 
2 3 

Due to the conservation of the total energy-momentum tensor the functions p{r) and 
s(r) are related to each other by the differential equation 

(5.9) 2 a£(r) ^ 2£(!l + Ml = . 

3 9r r dr 

Hereby p{r) describes the radial distribution of the 'pressure' inside the hadron, while 
s(r) is related to the distribution of the 'shear forces' [200]. 

Another important property which can be directly derived from the conservation of the 
energy-momentum tensor is the so-called stability condition. Integrating J d^r r''(Vi9'^ ) = 




by parts one finds that the pressure p{r) must satisfy the relation 



(5.10) 



dr r'^p{r) = 



The energy density Airr"^ pE{r)/MN = 4:Trr'^QQo{r)/MN normalized according to (5.7a) 
is shown in Fig. 11a as a function of r in the chiral limit of a vanishing pion mass and 
for the physical situation with a pion mass of 140 MeV. In this latter case the nucleon 
mass in the model is about 1250 MeV. This overestimate of the physical nucleon mass of 
about 300 MeV is typical for the soliton approach and its origin is well understood [201]. 
In the center of the nucleon the energy density is pe{0) — 1-70 GeV/fm'^ or 3.0 x 10^^ g 
cm~^, corresponding roughly to 13 times the equilibrium density of nuclear matter. As 
the pion mass decreases, the energy density is spread more widely. According to the role 
of the pion field in the xQSM, where one can associate the contribution of the discrete 
level to the quark core and the contribution of the negative continuum states to a pion 
cloud, this means that the range of the pion cloud increases and the nucleon becomes 
larger. Actually, the mean square radius (r|;) increases from 0.67 fm^ in the case of the 
physical pion to 0.79 fm^ in the chiral limit. With increasing pion mass up to 1.2 GeV 
this trend is confirmed with the nucleon becoming smaller and smaller [189]. 

The angular momentum density pj{r) is related to the 0ofe components of the static 
energy- momentum tensor as pj = e'^-'^ SiXjQi^k- The radial distribution of angular mo- 
mentum 4:Trr'^pj{r)/ Jn, normalized according to (5.7&) with the nucleon spin Jjv = ^, 
is shown in Fig. lib as a function of r for to^t — and 140 MeV. For any at small 
r one finds pj{r) oc . The mean square radius (rj) decreases with increasing m.^ [189] 
in agreement with the idea of a shrinking pion cloud. For a physical pion one finds 
(rj) = 1.32 fm^. At large r in the chiral limit pj{r) cx l/r"' such that (rj) diverges. 

Fig. 12a shows the pressure p{r) as function of r. In the physical situation p{r) takes 
its global maximum at r = with p(0) = 0.23GeV/fm^ = 3.7 • lO^-^Pa. This is 0(10- 
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Fig. 12. - (a) The pressure p(r) from the xQSM as function of r for mjr = and 140 MeV. (b) 
r'^p(r) as function of r at the physical value of (taken from [44]). 



100) higher than the pressure inside a neutron star. Then p(r) decreases monotonically 
(becoming zero at ro = 0.57 fm) till reaching its global minimum at rp_ = 0.72 fm, 
after which it increases monotonically remaining, however, always negative. The positive 
sign of the pressure for r < tq corresponds to the repulsion among quarks imposed by 
Pauli principle, while the negative sign in the region r > means attraction in agreement 
with the idea of a pion cloud responsible for binding the quarks to form the nucleon. The 
subtle balance between repulsion and attraction, ultimately producing a stable soliton, 
can be better appreciated from Fig. 12b where r^p(r) is shown as a function of r. The 
shaded regions have the same surface areas but opposite sign and cancel each other within 
numerical accuracy, thus satisfying the stability condition (5.10). 

Fig. 13 shows the distribution of the shear forces s(r) obtained by solving the differen- 
tial equation (5.9). The distribution of shear forces is always positive. For = 140 MeV 
it reaches a global maximum at r = 0.40 fm. The position of the maximum is weakly 
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Fig. 13. - The function s(r) describing the shear forces in the nucleon within the xQSM as a 
function of r (taken from [44]). 
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dependent on rriT^. At small r, sir) cx r^. 

5'3. Spin densities. - When ^ = and a; > 0, according to Refs. [62, 202] one defines 
three-dimensional densities 



(5.11) p{x,h,\,\N) = \ 



n{x, b^) + V>e^'S' — £'{x, b^) + AAjvW(a;, 6^) 



(5.12) p{x,b,s,S) 



n{x,b^)+s'S' Ht(x,62)- 



+ ^ (^S'S'ix, 6^) + [S!r{x, b^) + 2n'j.{x, 6^)]) 



4M2 



AbnT{x,b'') 



Mn 



+ s\2bV -b%)S^j^H'Ux,b') 



The distributions H, £, H, Ht, etc. are the Fourier transform of the corresponding GPDs 
as discussed in sect. 3.2 (see, e.g., Eq. (3.4)). In Eqs. (5.11) and (5.12) the shorthand 
notations 



(5.13) 



f' = — f 



d 



d 



have been used, and the two-dimensional antisymmetric tensor e*-' has been introduced 
with = — e^i — \ ^^id e^^ — e^^ = q. Roman indices are to be summed over. 

Since A+ = 0, p{x, b, A, Xn) and p{x, b, s, S) represent the probability to find a quark 
with longitudinal momentum fraction x and transverse position b either with light-cone 
helicity A (= ±1) in the nucleon with longitudinal polarization Ajv (= ±1) or with 
transverse spin s in the nucleon with transverse spin S. 

In Eq. (5.11) the first term with H describes the density of unpolarized quarks in the 
unpolarized proton. The term with £' introduces a sideways shift in such a density when 
the proton is transversely polarized, and the term with H reflects the difference in the 
density of quarks with helicity equal or opposite to the proton helicity. 

In the three lines of Eq. (5.12) one may distinguish the three contributions correspond- 
ing to monopolc. dipolc and quadrupole structures. The unpolarized quark density ^Ti. 
in the monopole structure is modified by the chiral-odd terms with Ht and A^Ht when 
both the quark and the proton are transversely polarized. Responsible for the dipole 
structure is either the same chiral-even contribution with £' from the transversely polar- 
ized proton appearing in the longitudinal spin distribution (5.11) or the chiral-odd con- 
tribution with £^ + 2H'rp from the transversely polarized quarks or both. The quadrupole 
term with H'^ is present only when both quark and proton are transversely polarized. 

Lattice calculations accessing the lowest two a;-moments of the transverse spin den- 
sities of quarks in the nucleon have recently been presented [99], and impact parameter 
dependent parton distributions in phenomenological models of hadron LCWFs have been 
studied in Ref. [203]. 

In Ref. [204] the first a;-moments of the spin distributions 



(5.14) p{b,X,XN) 



J dxp{x, 



6, A, Ajv), p{b,s,S)= / dxp{x,b,s,S) 
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Fig. 14. - The monopole contribution ^Ti (left) for unpolarized quarks, the dipole contribution 
— \sxhy{£'j- + 2'H'j^) / Mn (middle) for (transversely) x-polarized quarks, and the sum of both 
(right) in an unpolarized proton. The upper (lower) row gives the results for up (down) quarks 
(taken from [204]). 



have boon studied as functions of the transverse position and different quark and proton 
polarizations taking advantage of the overlap representation of LCWFs that was origi- 
nally proposed in Refs. [57, 141] and successfully applied to GPDs [88, 89, 91, 205, 146]. 

With such a model only valence quarks are considered. Therefore, the integrals in 
Eq. (5.14) are restricted to < a; < 1. In any case, for the unpolarized density this is not 
a limitation as the sea quark contribution drops out of the first moment of the GPD H. 
For the spin-dependent densities this restriction is not dramatic because valence quarks 
are known to dominate at large and intermediate x {x> 0.2). 

In Figs. 14 and 15 the distorting effect of the dipole terms due to the transverse 
spin distributions on the monopole terms corresponding to spin densities for unpolarized 
quarks in an unpolarized target is shown. If one multiplies the up and down monopole 
terms by the quark charge eq and sums over quark flavours, one obtains the nucleon 
parton charge density in transverse space, as discussed in sect. 5'1. 

For transversely polarized quarks in an unpolarized proton the dipole contribution 
introduces a large distortion perpendicular to both the quark spin and the momentum of 
the proton (Fig. 14). Evidently, quarks in this situation also have a transverse component 
of orbital angular momentum. This effect has been related [40, 206] to a nonvanishing 
Boer-Mulders function [149] which describes the correlation between intrinsic trans- 
verse momentum and transverse spin of quarks. Such a distortion reflects the large value 
of the anomalous tensor magnetic moment kt for both flavours, i.e. = 3.98 and 

= 2.60, to be compared with the values ~ 3.0 and ^f, « 1.9 of Ref. [99] due to a 
positive combination £t + 2Ht- Since kt ~ —hj-, the results of Ref. [204] confirm the 
conjecture of Refs. [40, 206] that hj- is large and negative both for up and down quarks. 

As also noticed in Refs. [13, 99] the large anomalous magnetic moments k"'** are 
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Fig. 15. - With unpolarized quarks the monopole contribution (left) for an unpolarized 
proton, the dipole contribution —^SxbyS' /AIn for a (transversely) i-polarized proton, and the 
sum of both (right). The upper (lower) row gives the results for up (down) quarks (taken 
from [204]). 
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Fig. 16. - Densities of unpolarized quarks in a transversely polarized nucleon (left) and trans- 
versely polarized quarks in an unpolarized nucleon (right) for up (upper plots) and down 
(lower plots) quarks. The nucleon and quark spins are oriented in the transverse plane as 
indicated, where the inner arrow represents the quark and the outer arrow the nucleon spin 
(taken from [99]). 
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Fig. 17. - The monopole contribution ^SxSxCHt — A[,'Wt/4M^) (left) and the quadrupole 
contribution ^SxSx{bl — by)'H''^/M% (right) for i-polarized quarks in a nucleon also polarized 
along X. The upper (lower) row gives the results for up (down) quarks (taken from [204]). 



responsible for the dipole distortion produced in the case of unpolarized quarks in 
transversely polarized nucleons (Fig. 15). With the model of Ref. [204], k" = 1.86 
and k'^ = —1.57, to be compared with the vahics ~ 1.673 and n'^ = —2.033 de- 
rived from data. This effect can serve as a dynamical explanation of a nonvanishing 
Sivers function [63] which measures the correlation between the intrinsic quark 
transverse momentum and the transverse nucleon spin. Such results, with the oppo- 
site shift of up and down quark spin distributions imply an opposite sign of /fi^ for 
up and down quarks [61, 64] as confirmed by the recent observation of the HERMES 
collaboration [207]. 

The results in Figs. 14 and 15 are in qualitative agreement with those obtained in 
lattice calculations [99] , where strongly distorted spin densities for transversely polarized 
quarks in an unpolarized nucleon have been found (Fig. 16). One observes that the 
sideways distortion for down quarks is about twice as strong as for up quarks, even if the 
anomalous magnetic moment and the anomalous tensor magnetic moment have 
about the same magnitude. This is because the monopole distribution for up quarks is 
twice as large as for down quarks, and therefore adding the dipole contribution results 
in a larger distortion for down quarks than for up quarks. 

In Fig. 17 quarks and proton all transversely polarized along x are considered. Quite 
remarkably, opposite signs of the up and down quadrupole and spin dependent monopole 
terms are found as a consequence of the opposite sign for up and down quarks of the 
X dependence of Ht and Ti.T predicted by the model [91]. The quadrupole distribution 
for up quark is more spread than for down quark, but it is much smaller, resulting in an 
average distortion equal to —0.04 to be compared with the value 0.07 for down quark. 
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Fig. 18. - The total spin distribution as a sum of monopole, dipolc and quadrupole terms, for 
x-polarized quarks in a proton also polarized along x; left (right) panel for up (down) quarks 
(taken from [204]). 



The total spin distribution for quarks and proton transversely polarized along x is 
shown in Fig. 18 as the result of summing for each flavour the two monopole contribu- 
tions in the left panels of Figs. 14 and 17, the two dipolc contributions on the middle 
panels of Figs. 14 and 15 and the quadrupole contribution of the right panel in Fig. 17. 
For up quarks this distribution is almost axially symmetrical around a position slightly 
shifted in the y-direction (Fig. 18). This is a consequence of the dominating role of 
the monopole terms, the compensating action of the two dipole terms and the small 
quadrupole contribution. In contrast the down quark spin distribution has a much lower 




A (to) 




-I'jr OA OS 



Km 







>0.J4 




t 014 




■ 0.2a 




J 0.22 


E 








m 




□ 


<-OLU 


□ 


<-lll.U 



■ •:11.0s 

[ I t-rMO 
I I <-(l<OS 



Fig. 19. - The dipole contribution ^SybxS' (left), the total dipolc contribution ^[SybxS' — 
Sxhy[£T + 2'H't)/Mn] (middle) and the quadrupole contribution SxSyhxhyliT / M% (right) for 
a-polarized quarks in a nuclcon transversely polarized in the y-direction. The upper (lower) row 
gives the results for up (down) quarks (taken from [204]). 
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Fig. 20. - The total spin distribution, as a sum of monopole, dipole and quadrupole terms, for 
a-polarized quarks in a proton transversely polarized in the y-direction; left (right) panel for up 
(down) quarks (taken from [204]). 



size and shows a strong and symmetric deformation about the y-axis stretching along 
the same x-direction of the quark and proton polarization. 

The resulting transverse shift of both up and down distributions in Fig. 18 is suggest- 
ing the presence of an effective transverse quark orbital angular momentum introduced in 
the LCWFs by the Melosh rotation required to transform the canonical spin to the light- 
front spin. Due to the shift in the positive y-direction transverse quark spin and trans- 




Fig. 21. - The monopole contribution ^TC (left) and ^Ti. (right), and their sum corresponding 
to the spin distribution for quark polarized in the longitudinal direction, parallel to the proton 
helicity. The upper (lower) row gives the results for up (down) quarks (taken from [204]). 
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verse quark orbital angular momentum seem to be aligned along the same a;-direction of 
the proton polarization. 

In Figs. 19 and 20 results are given for .x-polarizcd quarks in a proton polarized along y. 
The distortion due to the dipole contribution ^Syb^S' in Fig. 19 is rotated with respect to 
the case shown in Fig. 14, but the origin of opposite shift for up and down quarks is always 
the opposite sign of the anomalous magnetic moments k""''. The total dipole distortion 
in Fig. 19 is obtained by considering also the second dipole term —^Sxby{£'j, + 2'H'j,)/MM 
displayed in Fig. 14. The result is quite sizable, while the quadrupole term is rather 
small. Therefore, the total resulting distortion of the spin density (Fig. 20) is due to the 
dipole terms with a small contribution from the quadrupole terms. Correspondingly, the 
quark orbital angular momentum has positive x and y components for up quarks, and 
positive X and negative y components for down quarks. Here as well as in Fig. 18 the 
quark orbital angular momentum is entirely generated by the Melosh rotations. 

Finally, the case of quark polarization parallel to the proton helicity is considered in 
Fig. 21. Here only monopole terms occur (see Eq. (5.11)) and their role was first discussed 
in Ref. [62]. The opposite sign of ^Ti for up and down quarks is responsible for quite a 
different radial distribution of the axially symmetric spin density. Since in the forward 
limit the GPD Ti. reduces to the helicity distribution lS.q{x), H{x, ^ — 0,t — 0) — Aq{x), 
this difference ultimately reflects the opposite behaviour of the helicity distributions and 
the opposite sign of the axial-vector coupling constants g^'^ of up and down quarks (see 
also Ref. [89]). 

6. — Status of experimental investigation 

Experimentally, the GPDs can be accessed in exclusive measurements such as hard 
exclusive meson leptoproduction eA'' e'MN' [208, 12, 115], clcctroproduction of the 
photon eA^ e'N'j which is sensitive to the DVCS amplitude [8, 10], photoproduction 
of a lepton pair 'jN — > ii'N' or time-like Compton scattering [209] and electroproduction 
of a lepton pair eN e'N'fl' or double deeply virtual Compton scattering [210, 211]. 
At leading order the DVCS and meson production amplitudes can be described in terms 
of the so-called handbag diagram (Fig. 22), where the lower blob represents the involved 
GPDs. In the asymptotic limit of the corresponding differential cross sections scale 
as for DVCS and for meson production. These dependences are strong 
experimental signatures that the appropriate regime is reached in order to rely on the 
validity of the handbag description and the interpretation of data in terms of GPDs. 

The GPDs depend on three variables: x, £, and t. However, only two of them are 
accessible experimentally in lepton scattering, i.e. ^, fully defined by detecting the scat- 
tered lepton, and t, fully defined by detecting either the recoil proton or the outgoing 
photon or meson. Due to the loop in the handbag diagram (Fig. 22), x is integrated over 
and the GPDs enter the integral with a weighting function given by the propagator of 
the quark between the incoming virtual photon and the outgoing photon or meson. This 
means that the amplitude of the process is proportional to the convolution 

J -I (,-x + ie J_i ^-a; 

where F is any of the GPDs. The imaginary part of the amplitude is thus given directly 
by GPDs at the special point a; = ^, whereas the real part is sensitive to all x. However, 
also in the latter case the most relevant region of x is determined by ^. Only when 
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Fig. 22. - Handbag diagrams: (a) for DVCS and (b) for exclusive meson production. 



an observable is proportional to the imaginary part of the amplitude, like in the case 
for instance of the beam-spin asymmetry in DVCS, one actually measures directly the 
GPDs at some specific value .t = ^. It is therefore a non-trivial task to extract GPDs 
from data by a dcconvolution procedure, so that one has to rely on models and fitting 
procedures in a global study of several observables [212]. In addition, the effects of other 
mechanisms, such as Regge exchange contributions particularly important at large 
and small t [213, 214, 215], have to be considered. 

Hard exclusive meson production is harder to describe quantitatively than other pro- 
cesses, but it provides important complementary information such as gluon distributions 
that enter the Compton amplitude only at NLO order in the strong coupling constant a^. 
On the other hand, Compton scattering, especially double DVCS which provides a more 
direct means to measure quark GPDs and their dependence on the different kinematical 
variables, is much more demanding from an experimental point of view. 

Exploratory and first dedicated experiments have been performed in recent years at 
the HERA collider (27 GeV electrons/positrons against 800-900 GeV protons) with the 
HI and ZEUS collaborations, at DESY with the fixed-target HERMES detector and the 
27 GeV electron/positron beam, at JLab in Hall A and B with the 6 GeV electron beam. 
Other measurements are in program, as e.g. at CERN with the COMPASS experiment 
working with the high-energy muon beam. In the following a brief review of the available 
data is presented. 

6"1. Hard exclusive meson production. - Exclusive leptoproduction of light vector 

mesons, j*p Vp {V = p^,uj,(f)), at low virtuality of the virtual photon and for 
7*p cm. energies W ^ GeV exhibits features typical of soft diffraction processes: a 
weak dependence on W and a differential cross section that falls exponentially with —t 
at low jtj values (see, e.g., [216, 217, 218, 219] and references therein). These features are 
consistent with the expectations of the vector-meson dominance model and the colour- 
dipole model according to which the photon fluctuates into a vector meson which then 
scatters elastically from the proton. In contrast, vector-meson production at high Q"^ 
(hard diffractive leptoproduction) was proposed as an important source of information 
of the gluon dynamics [220, 11]. The production mechanism was assumed to proceed in 
three steps: the virtual photon fluctuates into a qq state; the qq pair scatters on the proton 
target via exchange of a gluon pair; and the scattered qq pair turns into a vector meson 
(Fig. 23a). It was also observed that longitudinally polarized photons with sufficiently 
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Fig. 23. - Gluon (a) and quark (b) GPDs contributing to the amplitude for meson production. 



large produce small-size qq pairs configurations [208]. This made possible to apply 
perturbativc QCD to calculate the cross section when the vector meson is produced in 
the longitudinally polarized state by longitudinally polarized photons at small x. In this 
case a factorization theorem for the process amplitude can be established [208]. The 
factorization theorem was proved in the coUincar approximation for the full range of x 
and the general case of hard exclusive electroproduction of any meson [12] . Contributions 
from transverse photons are suppressed by at least one power of 1 /Q with respect to the 
longitudinal ones [12, 221, 222]. The analysis of the cross section for meson production 
with target polarization can be found in Refs. [223, 224]. 

In the scaling limit, — > oo, at fixed Bjorken variable xb = Q'^/'2{p ■ q) and fixed 
invariant momentum transfer t to the proton, the amplitude of vector- meson production is 
given by a convolution of the soft parts of the process, consisting of both gluon and quark 
GPDs [11, 12] and the nonperturbative qq distribution amplitude for meson formation, 
with the perturbativcly calculable hard-scattering kernels (Fig. 23). The kernels are 
known to NLO, i.e. to order [225]. The description is restricted to sufficiently large 

but can be used for both small and large xb, thus providing a common framework 
for analyzing both collider and fixed-target data [226]. In the coUinear factorization 
framework NLO corrections are not negligible in the xb range relevant for experiments 
at HERMES, JLab and COMPASS [227]. 

Hard exclusive meson productions are rather complex to analyze as they contain 
nonperturbative information on both the target and the produced meson. Nevertheless 
they offer the possibility to disentangle different GPDs. For mesons with natural parity 
P = (— )'^, such as vector mesons, both quark and gluon GPDs in general contribute 
at leading order in ag, whereas for mesons with unnatural parity P = (— )'^^^, such as 
pseudoscalar mesons, only quark contributions appear. The quantum numbers of the 
produced meson select different flavour combinations of quark GPDs. At leading twist 
exclusive vector- meson production is only sensitive to unpolarized GPDs (H'^ and i?''), 
whereas polarized GPDs (7?' and E'^) are involved in pseudoscalar meson production 
without the need for a polarized target or beam [228] . The gluon GPDs and can 
be accessed in neutral vector- meson production [115, 116, 229, 230, 231]. 

Chirally odd GPDs cannot be accessed in exclusive electroproduction of a single vector 
meson [221, 222, 232]. They can be accessed in diffractive electroproduction of two vector 
mesons [233]. 

First results on hard diffractive p° production were obtained in the collider experi- 
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ments HI [234] and ZEUS [235, 236, 226] by looking at the W, Q"^ , and t dependence 
of the cross section in terms of the invariant mass distribution of the two decay pions. 
Evidence for a dominant longitudinal contribution, increasing with Q^, was found [235]. 
At large the cross section also develops a stronger W dependence than that expected 
from the behaviour of elastic and total hadron-hadron cross sections [234] . Therefore, no 
comparison with models based on GPDs was possible. Also the recent high-precision data 
of the ZEUS collaboration [226] have been compared to various theoretical predictions, 
none of which are able to reproduce all the features of the data. 

With fixed hydrogen target [237, 238] the meson decay into tt+tt" was used to 
identify the reaction of interest. The longitudinal part of the cross section was extracted 
from data [237] or separated under the assumption of s-channel helicity conservation, 
i.e. that transitions with the same helicity for photon and meson are larger than those 
changing the helicity [238] . A fair description of data is obtained with the phenomenolog- 
ical GPDs parametrization of Refs. [115, 116] with the dominating quark contribution, 
while the gluon contribution starts to contribute significantly only at = 4 GeV^ and 
W >IQ GeV. 

Hard exclusive electroproduction of tt+tt" pairs has also been measured [239] . In this 

case both two-gluon and quark-antiquark exchange mechanisms contribute. The resonant 
7r+7r~ pair production via longitudinal p° (isospin / = 1, total angular momentum J = 
1,3,... and C-parity C = —1) or /-meson {I — 0, J = 0,2, . . ., C — +1) decay is in 
competition with the two-gluon channel that gives rise to pion pairs with the quantum 
numbers of the p-meson family only. Therefore a sizable admixture of isoscalar and 
isovector pion pairs is obtained. 

The cross section of hard exclusive electroproduction of tt"*" mesons was measured 
by the HERMES collaboration as a function of t and [240, 241, 242]. A model 
calculation based on GPDs with power corrections [116] is in fair agreement with the 
data at low values of \t\. However, in the experiment separation of the longitudinal and 
transverse photon contributions to the cross section was not achieved. The t and 
dependence of the cross section was nicely described by a model calculation based on the 
Regge formalism [213]. This good agreement is mainly due to the dominant pion-pole 
contribution which enters also in the GPD E. 

The exclusive oj electroproduction off the proton was studied at JLab [243] in a large 
kinematical domain above the nucleon resonance region. The f-channel 7r° exchange, 
which is mostly due to transverse photons, seems to dominate the reaction even for 
as large as 5 GeV^. Thus the GPDs formalism is hardly applicable in this case. 

6'2. DVCS. - After some seminal work [8, 10, 9] a factorization theorem was proved 

valid also for DVCS to all orders in perturbation theory [71, 244]. Compared to exclusive 
meson production DVCS is simpler to describe because the composite meson in the final 
state is replaced by an elementary particle, the photon, and thus there is no meson wave 
function in the factorization formula. Perturbative contributions to DVCS have been 
worked out to NLO accuracy [245, 246, 247]. A complete analysis in the twist-three 
approximation of the theory of DVCS was developed in [248] , where an exhaustive set 
of analytical results was given for the cross section for all possible hadron and lepton 
polarizations involved. A detailed discussion of lepton scattering on longitudinally or 
transversely polarized protons is given in [223]. 

The reaction ep — > e'p'j receives contributions from the DVCS process and the purely 
electromagnetic Bethe-Heitler (BH) process where the photon is emitted from the initial 
or the final electron (Fig. 24). While the BH cross section in most of the accessible 
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Fig. 24. - (a) Compton scattering; (b) and (c) Bethe-Heitler contributions. 



kinematic regions is much larger than the DVCS cross section, the DVCS contribution can 
be measured through interference of DVCS and BH amphtudes. According to theory [249, 
248], the leading-order and leading-twist DVCS-BH interference is proportional to the 
sum of two terms, one helicity independent {da) and one helicity dependent {dS), i.e. 



where ) denotes a negatively (positively) charged lepton with polarization P;, is the 
azimuthal angle of the produced photon with respect to the lepton scattering plane, and 
e is the polarization parameter of the virtual photon. The amplitude M.^'^ is given by a 
linear combination of the Dirac and Pauli form factors Fi and F2 together with Compton 
form factors H, H and £ that are convolutions of the GPDs H, H and E, respectively. 
dT, and da provide complementary information: dT, measures the imaginary part of 
providing direct access to GPDs at x = ^, whereas da determines the real part of the 
DVCS-BH interference and depends on the integral of GPDs over the full range of x. 
The beam-spin asymmetry, as the ratio of the difference to the sum of cross sections 
with a polarized lepton beam, is associated with the imaginary part of the DVCS-BH 
interference with a sin cj) modulation. The beam-charge asymmetry, as the ratio of the 
difference to the sum of cross sections with opposite charge of the incoming lepton beam, 
accesses the real part with a cos (j) modulation. 

First DVCS data on the proton have been reported by measuring the cross sec- 
tion in collider experiments at high energy by the HI [250, 251] and ZEUS collabora- 
tions [252, 253] and by looking at the beam-spin asymmetry in fixed-target experiments 
at lower energy with polarized lepton beams and unpolarized targets by the HERMES 
collaboration [254] and at JLab [255]. More precise DVCS data are now available on 
beam-spin asymmetry [256] as well as on the longitudinal target-spin asymmetry [257] 
from JLab and on transverse target-spin asymmetry from HERMES [258]. The lepton 
beam-charge asymmetry has also been measured at HERMES [259, 260] and HI [261], and 
will be a dedicated experiment in a future program at COMPASS with the high-energy 
muon beam at CERN [262] . Quite recently the neutron contribution to the DVCS off a 
deuterium target has been extracted from the helicity-dependent cross section measured 
at JLab in Hah A [263]. 

The consistent HI [250, 251] and ZEUS [252, 253] cross section results for up 
to 100 GeV^ and W up to 140 GeV are compatible with NLO calculations using GPDs 
parametrizations [264] as well as with a description in terms of colour dipole models [265, 
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266] which have been successful in describing both the inclusive and diffractive DIS cross 
sections at high energy. 

First results on the beam-spin asymmetry with the positron beam at HERMES [254] 
and the electron beam at JLab [255] confirmed a sincf) modulation in agreement with 
calculations based on GPDs parametrizations predicting a dominant contribution of the 
GPD H [115, 117]. 

The full potential of the beam-spin asymmetry has been first explored at JLab [256] 
by separating dS and da. The absence of dependence in the range 1.5 < < 2.3 
GeV^ of the sin cf) term of dE supports the twist-two dominance of the DVCS amplitude 
driven by the GPD H. The dT, data are in qualitative agreement with the predictions 
from the model of Refs. [116, 16, 74], but da is significantly underestimated. 

First measurements of the beam-charge asymmetry have been reported by the HER- 
MES collaboration using the electron and positron beams at HERA [259, 260]. A coscj) 
dependence has been observed in the range 0.03 < —t < 0.12 and xb — 0.1. The data are 
in agreement with the dual parametrization of GPDs [73] and the dominance of the GPD 
H with a large negative D-tcrms [16, 117]. However, the data allow sufficient freedom 
in modeling the unknown small- a; behaviour of the double distribution part of GPDs, so 
that it is also possible to describe them without the D-tevm [248] . 

The longitudinal target-spin asymmetry Aul in ep^ e'p'7, as the ratio of the differ- 
ence to the sum of cross sections with target polarization antiparallel and parallel to the 
beam direction, has been first measured at JLab [257] . A dominating sin (j) modulation 
with a large contribution from GPD H is observed, consistent with predictions based on 
the GPD formalism. 

The transverse target-spin asymmetry Aut associated with DVCS on the proton is 
measurable as the ratio of the difference to the sum of cross sections using an unpolarized 
lepton beam and a transversely polarized target proton at angles ^5 and (ps + n with 
respect to the lepton scattering plane [258]. The interest of such an experiment lies 
in the sensitivity of Ajjt to the GPD E entering the Compton form factor driving a 
cos^ modulation. In turn, E can be modeled using the total angular momentum 
carried by quarks in the nucleon [16]. Since the contributions of the up and down quarks 
are proportional to the corresponding squared charge, the down quark contribution is 
suppressed and Aut can give constraints on J" [113]. Model-dependent constraints can 
indeed be derived from experiment [258] in good agreement with lattice simulations [93, 
48]. Using the model of Ref. [16] the result is J"-F J''/2.9 = 0.42 ±0.21 ±0.06. The quite 
recent analysis of the beam-spin asymmetry on the deuteron [263] provides a correlated 
constraint on J" and J'^ from a fit to the neutron data {J'^ + J"/0.5 = 0.18 ± 0.14) 
intersecting the constraint from the proton in close vicinity of the lattice results. However, 
one should remind that the lattice results were obtained by neglecting disconnected 
diagram contributions which only drop out of the isovector combination J" — J'^ but can 
be relevant for the isoscalar combination J'^ + J'^. 

At the present stage of experimentation it is clear that DVCS is a quite selective pro- 
cess to study GPDs. However it is also clear that, in absence of any model-independent 
deconvolution procedure, one has to rely on some global fitting of data [77]. Therefore, 
a full experimental program aiming at the extraction of the individual GPDs requires to 
study several observables in different channels with dedicated experiments. 
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7. — Conclusions 

The study of the internal structure of the nucleon through the newly developed con- 
cept of GPDs has seen substantial progresses in the last decade. The physical content of 
the GPDs has fully been disclosed through the theoretical study of these functions both 
in momentum space, where they allow to study the momentum correlations of partons 
in the nucleon, and in impact parameter space, where they enable to map the spatial 
densities of the nucleon. At the same time, the rich spin structure of the GPDs provides 
ways to study aspects of the nucleon spin which are otherwise difficult to obtain, like for 
example the role of the orbital angular momentum of nucleon's constituents in building 
up the total spin. 

The challenging task of the theory is to develop models able to capture the dynamics 
responsible for these many non trivial features of hadron structure. Efforts have been 
made in several directions, some of them with more emphasis on the phenomenological 
implications, others within the framework of effective models aimed to give complemen- 
tary descriptions of the nucleon dynamics. 

Among the phenomenological approaches, we discussed different strategies which in- 
corporate the constraints of the parton distributions in the forward limit and of the 
electroweak form factors for the first moments in x of the GPDs. As a very promising 
tool to handle the phenomenology of the GPDs we also considered recent fitting proce- 
dures which use additional inputs from lattice QCD simulations for the GPD moments. 

As far as effective models are concerned, we reviewed the basic ideas and results 
obtained in the chiral quark-soliton model, light-front constituent-quark and meson- 
cloud models, with the aim to unravel the different role of the valence and sea quarks in 
the nucleon's structure. These hadronic models refer to low energy parametrizations of 
the GPDs and need to be supplemented by QCD evolution before the comparison with 
experimental results at larger scales. 

Ah initio calculations of the GPDs can be obtained through lattice QCD simulations. 
Lattice QCD methods have made steady progress during the last years and first results 
became recently available for moments of parton distributions, form factors and GPDs. 
Although the actual simulations refer to unrealistic large quark masses and only to the 
fiavour non-singlet quark contribution, they have provided important insights into the 
transverse structure of unpolarized nucleons, the lowest moments of polarized and tensor 
GPDs, and transverse spin densities of quarks in the nucleon. Significant investments 
are being made to improve the accuracy of the lattice methods and to extend the calcu- 
lations also to the gluon and the singlet-quark GPDs, holding the promise to provide a 
complementary tool to the experimental investigation of the GPDs. 

Experimental information about GPDs come from measurements of a variety of hard 
exclusive processes in eN scattering, which rests on the possibility to apply factoriza- 
tion theorems to unambiguously separate the short-distance information specific to the 
probe from the long-distance information contained in the GPDs, representing universal, 
process-independent quantities. 

Deeply virtual Compton scattering and meson production in eN scattering have been 
measured in several experiments at fixed-target facilities (HERMES at DESY, Jeffer- 
son Lab in Hall A and B) and the HERA collider. The most problematic task in such 
measurements is to disentangle the dependence of GPDs on the three variables x, ^ and 
t. This is becoming possible with the refinement of procedures which combine theo- 
retical studies, phenomenological parametrizations, non-perturbative calculations such 
as lattice QCD, and increasingly more accurate experimental data. New measurements 
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both for DVCS and meson production are in program at CERN with the COMPASS 
experiment working with the high-energy muon beam. The HI and ZEUS cohaborations 
will continue to analyze data on DVCS and vcctor-mcson production at high and 
small X. The JLab 12 GeV upgrade program will double the energy beam of the CEBAF 
continuous-beam electron scattering accelerator from 6 to 12 GeV. Thanks also to the 
upgrade in the detection equipment of the experimental halls [267], high quality data arc 
expected for exclusive processes in a larger Q^, x and t kinematical range. The de- 
pendence will provide stringent tests of the factorization theorem, while the possibility to 
scan the t dependence of the obscrvablcs as a function of x will allow to study the trans- 
verse profile of the nucleon as a function of quark light-cone momentum fraction [268]. 
To separate the different spin components of the GPDs, measurements of a variety of 
polarization obscrvablcs (beam and target spin) are also planned. Moreover, the flavour 
decomposition of the GPDs will come both from measurements of DVCS with proton 
and deuterium targets, and from meson production experiments. 

Finally, the construction of a new Electron Ion Collider (EIC) is currently under dis- 
cussion for a location in the USA [269] . The capabilities provided by this machine would 
complement the information from fixed-target experiments, oflfering the opportunity to 
map the gluon and quark GPDs at very small x. 
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